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Appendix B 


SUMMARY OF FEAP MODIFICATIONS FOR IMPROVED EFFICIENCY 


The major changes made to the FEAP program are listed below: 


1. The entire program was recompiled using FTN (OPT=2). Several FORTRAN errors were 
found and corrected. 


2. The PUNCH, PUNCHB, TAPE8 and TAPE12 files were removed. Buffer sizes were 
increased to 3075 for TAPE7 and TAPE10 and to 2050 for TAPE13 and TAPE 14. 
TAPE11 becomes a direct access (QUICKIO) file with a buffer size of 65 words. QUICKIO 
isa FORTRAN callable subroutine which decreases I/O time. 


3. All files were binary blocked using FTNBIN. 


4. The global stiffness matrix is saved on a direct-access file in subroutine MAKEA. During 
forward reduction in USOL, no additional files are used for spill-over processing. All 
intermediate storage required is written in-place as an appropriate record on the QUICKIO 
file. This reduces the appropriate record on the QUICKIO file, which in turn reduces the 
amount of scratch storage required by more than half. 


5. USOL was completely rewritten to allow direct-access equation processing. A COMPASS 
routine was included to increase the solution efficiency. 


6. RESOLVE was rewritten to accommodate direct-access processing during the back- 
substitution phase. 


7. Subroutines PRTMAT, COMPACT, CKQUAD, PACKD, SLDO4, and QSHP8 were 
removed. 


8. Subroutines FELMT, MESHCK, TSOLVE, and RESET were modified to accommodate 
the node resequencing (bandwidth minimizing) changes. Primarily, this involved adding 
named COMMON/MINBW/NEW(1500) and internally altering nodal references. Note that 
fixing NEW at 1,500 words inherently limits the program to 1500 nodes maximum. This 
COMMON would have to be extended if larger problems are anticipated, 


9. Two input changes were made, neither of which affect the execution of existing decks; 
1.€., no actual changes were made to ir put processing unless the new features are desired. 
An optional first card may be included with the entry ISZDT. Blank COMMON will be 

set to this value allowing execution at reduced core size tor smaller problems. The entry 

is read using IS format. The default value is 58000, The second option is activated by 
entering “1” in ce 79 of the FEAP card. If activated, the program expects to encounter 
the PREFEAP punched output immediately tollowing the FEAP card. PREFEAP punched 
output contains the optionally reordered nodal point numbering system. 
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Appendix C 


USER INSTRUCTIONS FOR THE HIM CODE 


|} CONTROL CARDS 


Card 1. (15) 


Columns Variable 


3-3 ISZDT 


Card 2. (6X, 12A6) 


Columns Variable 


7-10 


Card Group 2a. If Card 2 contains a 1 in ce 79, the program expects to encounter a group 
of cards specifying an optimally reordered nodal point numbering system prior to encountering 


Card 3. 


Card 3. (15,1X,3A6) 


Columns Variable 
1-5 NOIM 
7-12 

13-18 

19-24 


Card 4. (15,1X,3A6) 


Columns Variable 
1-5 NDF 
7-12 

13-18 

19-24 


Description 


Blank COMMON will be set to this 
value. Default value is 58000 


Description 


Must contain the word FEAP 73 
Output page header 


Enter 1 if optimized nodal point 
numbering is employed. Otherwise 
leave blank 


Description 


Spatial dimension of problem. Use 3 
X coordinate designation 
Y coordinate designation 


Z coordinate designation 


Description 


Number of unknown displacements per 
node. Use 3. 


X displacement designation 
Y displacement designation 


Z displacement designation 


Note 


(1) 


Note 


(2) 
(3) 


Note 


Note 


Card 5. (15,10X,15,30X,3F 10.0) 


Columns Variable Description Note 
1-5 NEN Maximum number of nodes connected 


to any one element. Use 8. 


16-20 MBAN Maximum expected half-bandwidth (4) 
51-60 CON3 (1) User defined X scale factor (5) 
61-70 CON3 (2) User defined Y scale factor 
71-80 CON 3 (3) User defined Z scale factor 
4 
ra 


Notes 


(1) For typical HIM code simulations this card can be left blank thereby specifying 58000 | 
as the blank COMMON allocation. 


(2) Any description that the user wishes to have printed at the top of each output page 
can be entered. 


(3) Large, three-dimensional, dynamic problems can be expensive. A bandwidth minimizer 
has been employed to resequence the nodal point numbering scheme and optimally reorder 
the system of equations to provide economy in computational costs. A set of cards has been 
provided for this purpose which map the old and new numbering schemes. The mapping is 
transparent to the user; he refers always to the old system. Use of this option is highly 
recommended, 


(4) Typically this value is 150. The default value is 100 and is not sufficient. Therefore 
150 should be specified. 


(5) These scale factors can be used to alter the size and shape of the HIM skull configuration 
as provided. Their default values are unity. Also these scale factors can be employed to 
change the provided geometrical configurations from inches to any appropriate unit of 
measure; i.€., centimeters. 


I} MATERIAL PROPERTY SPECIFICATION 


Card 1. (IS, 1X,12A6) 


Columns Variable Description Note 
1-5 NUMMAT Number of different material tT) 


characterizations to follow. 
7-12 Must contain the word MATERI 


One group of cards follows for each material. There must be NUMMAT card groups. There 
are at least two cards for each group; card (a) and card (b). For viscoelastic characterization, 
additional cards are required within the group as described below. 


Card (a). (15,1X,15,11A6) 


Columns Variable Description Note 
1-5 Material number (1 to NUMMAT) 
7-11 . Element type (2) 
12-77 Alphanumeric description of each 

material 


Card (b). (215,2F 10.0) 


Columns Variable Description Note 
1-5 NG Number of viscoelastic shear terms. If (3) 


not viscoelastic, leave blank. 


6-10 NK Number of viscoelastic bulk terms. If 
not viscoelastic, leave blank. 


11-20 The infinite shear modulus or elastic 
shear modulus, Go. | 
21-30 The infinite bulk modulus or elastic 


bulk modulus, Ko. 


31-40 Mass density 


Card (c). (8F 10.) Omit this and the next card if material 
is not viscoelastic. 


Loe nas Webebic Deseeipeon Ste. 
1-10 First shear modulus term Gq 

11-20 First relaxation time term, A 4 

21-30 Second shear modulus term, G9 


etc.,until NG terms are identified, 


Card (d). (8F 10.) 


Columns Variable Description Note 
1-10 First bulk modulus term, Ky 

11-20 First bulk relaxation time term, By 

21-30 Second bulk modulus term, K 


etc., until NK terms are identified. 


ed 
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Notes 


(1) NUMMAT is usually assigned a value of 5. This ‘lows the option for a different material 
specification tor each of the following skull-brain regions. 


Material Number Region 
I Outer table bone 
7. Diploe 
3 Inner ‘able bone 


Subarachnoid space 
Brain 


an 


(2) For the elastic skull materials (1, 2, and 3) the word ELM21 should be entered. For 
the subarachnoid space material (4) and brain material (5), an element possessing reduced 
integration is recommended, and the word ELM22 should be entered. Otherwise ELM21 


may be entered for all materials. 


(3) Viscoelastic material characterization is specified by the following equations: 


NG 

Git) Go + > Se" 
i=] 

= B 

2 -t/B; 

K(t) = Ky + Kic 1 
J=1 


Recommended units are the inch-pound-second system, although any consistent system 


can be used. (See note 1(5)). 
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Hl SKULL-BRAIN GEOMETRICAL INPUT 


Card 1. (15,1X,A6) 


Columns Variable 
1-5 NUMNP 
7-12 


Description 


Number of nodal points in the skull- 
brain system 


Must contain the word NODAL 


There must follow NUMNP nodal point cards in numerical order. 


Typical Nodal Point Card (15,115,3F 10.0) 


Columns Variable 


Description 


Nodal point number 


Insert O for free node or 1 for restraint 
in X direction 


Insert O for free node or 1 for restraint 
in Y direction 


Insert O for free node or 1 for restraint 
in Z direction 


X coordinate value 
Y coordinate value 


Z coordinate value 


Note 


(1) 


Note 


(2) 


Notes 


(1) For the six-cubed, half-skull model configuration provided, this value is 716. This is 
the model which has been most useful in the development of the HIM code. A complete 
set of nodal point cards whose coordinates are entered in inches has been provided tor 
this configuration. 


(2) The boundary conditions to be used foi the HIM simulation are activated in ce 15, 16, 
and 17. Entering a 0 allows the node to move in any one of three possible directions. A 
zero in cc 15, 16, and 17 means the node ts completely free. A 1 means that the node's 
displacement is constrained to move in a manner to be specified in subsequent input cards 
(see Section VI). The sets of nodal point cards which have been provided do not have any 
boundary condition preset with the following exceptions. The nodal point cards contain « 

1 in ce 16 for all nodal points located in the midsagittal plane. Thus Y component displace 
ments for half skull simulation are automatically constrained to zero in the madsagittal plane 
for the user’s convenience. 
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[Vo ELEMENT CARDS 


Card 1, (15. 1X,A6) 


Columns Variable 
1-5 NUMEL 
7-12 


Description 


Number of elements in the skull-brain 


system 


Must contain the word ELEMEN 


There must follow 2 cards for each element, taking the elements in numerical order, 


Typical Element Card Set 


Card (a) (215) 


Columns Variable 
1-5 
6-10 


Card (b) (814) 


Columns Variable 


1-4 
5-8 


9-12 


29-32 


Description 
Element number 


Material number 


Description 
Node 1 
Node 2 
Node 3 


Node 8 


Note 


(2) 


Note 


seinen 


Notes 


(1) For the six-cubed, half-skull model configuration, the value is 555. A complete set of 
clement cards has been provided for this configuration, 


(2) A material number for which properties have been identified in Section If is entered. 
The corresponding element is then assigned those properties. 
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V INITIAL CONDITIONS 


These cards can be omitted if the initial displacement velocity and acceleration of the 
skull-brain system are zero 


Card 1. (15, 1X,A6) 


Columns Variable Jescription Note 


1-5 NICD Number of initial condition vectors. (1) 
Specity a “3 


7-12 Must conta: « the word INITIA 


Card 2. (6X,2A6) Repeat NICD times 


Columns Variable Description Note 
7-18 Descriptive title for initial (2) 
conditions 


Card 3. (215,7F 10.0) Repeat for each degree of freedom 


Columns Variable Description Note 
1-5 Degree of freedom number (3) 


for X; 3*NODE -2 = entry 
for Y; 3 NODE -1 = entry 
tor Z; 3*« NODE -O = entry 


6-10 Generator increment (4) 
11-20 Initial displacement value 

21-30 Initial velocity value 

31-40 Initial acceleration value 


Card 4. (15) 


6) serene ser ee ree. 


Columns Variable Description Note 
1-5 Largest node times 3 plus 1 (5) 
(NUMNP* 3+1) 
13 


ine ell - — —_ — a 


a - 


Notes 


(1) When initial conditions of the skull-brain system are specified, the initial displacement, 
velocity, and acceleration (three vectors) are supplied as input. 


(2) Three cards follow Card 1 and merely describe each of the three initial condition vectors. 
These descriptions will appear above each vector in the subsequent printout. 


(3) Any or all nodes may be assigned an initial condition in the X, Y, and Z directions. 
But the initial condition must correspond to a degree of freedom number associated with 
the node and not the node number itself. The indicated equation is employed to obtain 
the degree of freedom number. 


(4) The generator increment is used when the initial condition is uniformly applied over 

the entire skull-brain system, and this is generally the case. It avoids the necessity for inputting 
all but a few cards. For example, if the skull-brain system is initially moving at 240 in./sec 

in the +X direction, two cards are required to assign the initial velocity to all nodes. The 

first would be for node 1 and the second for node NUMNP (largest node number in the 
system). The first card would contain a degree of freedom number of 1(3* 1-2), a generator 
increment of 2 (assigns the initial condition to only the degree of freedom numbers associated 
with the X-direction), and a 240 in ce 21-30. The second card would contain the value of 

3* NUMNP-2 in ce 1-5 and 240 in ce 21-30. 


(5) This card ends the initial condition data. 


. 
i 


Sg SE | tre ree. 


VI PROPORTIONAL LOAD OR DISPLACEMENT SPECIFICATION 
USING ANALYTICAL FUNCTIONS 


This section should be skipped if using point-by-point input loading. Such loading is 
covered in Section VIII. But this section is applicable to simulated dynamic loads or 
displacements which rise and subside according to a single time function (see Section VIII 
for the time function). They may all be assigned different magnitudes, but otherwise at 
any given instant the loads are proportional to one another. 


Card 1. (15, 1X,A6) 


Columns Variable Description Note 


1-5 Last node to which a force or (1) 
displacement is to be specified. 
Must contain the word FORCE. 


The following group of cards is composed of one card for each node to which a specified 
force or displacement is assigned. 


Typical Card (15,5 X,6F 10.0) 


Columns Variable Description Note 
1-5 Node to which a force or displacement 


is specified 


11-20 X component of force or displacement (2) 
21-30 Y component of force or displacement 
31-40 Z component of force or displacement 

15 


Notes 


(1) Merely means the largest node number to which a force or displacement ts applied 
Forces and displacements must be applied at nodal points. 


(2) This value is the ratio between the simulated peak load component at this node and 
the largest peak simulated load considering all nodes. It may also be a direction cosine. 

Only ratios from 0 to 1 are assigned here. The magnitude of the resultant for the largest 
peak load is prescribed in Section VUE (see Note (9)). 


; | 
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VIL PLOT CARD 


This card is to be used only when diverting computed response data to an output 
tape. This action does not suppress the printing of output data. 


Card 1. (15,1X,A6) 


Columns Variable Description Note 


1-5 Number of time steps plus one for (1) 
which data is output 


| 
Nm 


Must contain the word PLOT3D 
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Note 


(1) Data will be placed on an output tape for subsequent postprocessing for each time 
step up to, but not including the value indicated, Usually data for all time steps is desired. 
Pherefore this entry should be one more than the number of time steps prescribed in the 
dynamic analysis. See the variable NTS in the next section, 
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Vill IMPLICIT TIME-INTEGRATION SPECIFICATION 
his card is required on all static or dynamic runs. 


Card 1. (15,1X,A6) 


Columns Variable Description Note 
1-5 NSEQ Number of time sequences, usually 1 (1) 
7-12 Must contain the word IMPLIC for 


dynamic problems. Must contain the 
word VISCOE for static problems. 


NSEQ sets of the following cards are required 


Card (a) (F10.0,815,2F 10.0) 


Columns Variable Description Note 

1-10 DT Time step size (2) 

11-15 NTS Number of time steps in the sequence(s) (3) 
plus 1 

16-20 INT Print interval (4) 

Z1-25 NNI First node printed, generally 1 (5S) 

26-30 NNE Last node printed, generally NUMNP 

31-35 NEI First element stress printed, 


generally 1 
36-40 NEF Last element stress printed, 
generally NUMEL 
41-45 NPROP Number of proportional loads (6) 
46-50 NFORC Largest node to which a force or (7) 
displacement has been specified in 
a point-by-point specification of 
input loading 
51-60 BETA Newmark integration parameter, (8) 
usually specify 0.25 
61-70 DEI Newmark integration parameter, 
usually specify 0 
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Card (b). (215,7F 10.0) One for each NPROP. Use only when specifying analytical 


function input, Tf using point-by-point input skip to card (c). 


( olumns Variable Description Note 
1-5 Proportional load type, either 1 or 2 (9) 
6-10 K Exponent for type 2 

11-20 TMIN Start ume for loading 

21-30 PMAX End time for loading 

31-40 AO 

41-50 Al 

51-60 A2 

61-70 A3 

71-80 A4 


Card group (c) (15,5X,6F 10.0). This group of cards is comprised of one card for each node 


to which an input value is assigned. Node NFORC is the last node from which a card is 
provided. The group represents a loading vector, and generally there should be one group 
for each time step in the dynamic analysis. Omit this card group it NPROP is non-zero. 


Columns Variable Description Note 


1-5 Node to which force or displacement 


is specified 


11-20 X component value of force or 
displacement 

21-30 Y component value of force or 
displacement 

31-40 Z component value of force 


Card (3) (AS) last data card 


Columns Variable Description Note 
1-5 Must contain the word STOP 
20 
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Notes 


(1) Little experience has been gained with values other than 1. Potentially, NSEQ can be 
used to specify a broader range of different loading functions. (Card (a) describes how 
different funetions can be specified with NSEQ > 1). For example, one analysis could 
consist of two sequences where the first is the onset of loading requiring a small time step 
and the second is the loading decay which may be slower and during which a larger, more 


economical, time step can be employed 


(2) This value should be chosen with caution. [It generally depends on the rise time and 
should be no greater than one-fifth of that value. A value of onc-tenth would be sufficiently 


small 


(3) This is the number of time steps plus 1 for which the dynamic analysis is to proceed 
Generally, it is chosen so that the entire rise and decay of the loading function can be 


simulated with the prescribed time step. 


(4) Several inches of paper may be required, depending upon NTS, for the printing of output 
information. A blank in ce 16-20 will cause the printing of data for each time step, issuing 


a 1 will cause the printing of data for every other time step, ete. 


(5) The next four entries also provide for reductions in the amount of output. Tf a few nodes 
and clements can be identified as key check points in the analysis it is only necessary to 
monitor their output. However, data elsewhere should at least be saved on tape (see Section 
VII). 


(6) This value is employed to distinguish among the specification of different input load 
functions in the dynamic analysis. Kor example if it is desired that the prescribed load vary 
as a haversine function and subsequently change to a linear function, then NPROP=2. This 
is also true of a triangular function. If a single function, such as a ramp function or haversine 
function, is desired then NPROP=1. [fa point-by-point input loading is desired then NPROP 


is lett blank. 


(7) This entry is for point-by-point load specification only. If proportioned load functions 
are employed, NFORC is left blank. 


(8) BETA and DEL are the Newmark integration parameters (implicit). Investigators have 
determined that BETA = 0.25 and DEL = 0 are good choices. If DEL > O artificial positive 


damping ts introduced and if DEL < 0 artificial negative damping is introduced. The 


numerical integration will be unconditionally stable if 


4 (DEL + 1)2 
: 16 


BETA 


The higher and least important nodes may be damped out by choosing BETA = 0.276 and 
DEL = 0.05. 


ONT § ORS RR 


(9) Proportional load type 1 isa polynomial function, 
fit) = AO + Alet + Ader? + A3at? + Adet® 
Proportional load type 2 isa tr inscendental function, 
f(t) AO#SinK(Atet) + A2e*cosK(A3et) + A4 
\ haversine function has been found to be particularly useful in the specification of impact 
forces, F(t) = AO«SinK(A Lt), where AO is the peak force magnitude, A1 is pi dis ided 


by the load duration in seconds, and K = 2. Force magnitude must be halved when employing 


a half-skull model. 


ee 
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Appendix D 


PLOW DIAGRAM FOR MAIN SUBROUTINES 
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Set size of blank common DI USZDT) 


PRAP72 


MESH(1) 


Main driver routine, calls major subroutines allocates storage 


MATERIL, sets up material properties 


NODAL, or, GENERA, sets up nodal data 


MESH (2) 


MESHI(3) 


ELEMEN, sets up clement data 


BLOCK, sets up nodal and element data 


BRKGEN 


FORCE, sets up nodal force vector 

.——_—_————_ BOUNDA, overwrite nodal boundary conditions 

——_————__ BLOADS, surface pressure loads 

—_ —— _ FLOADS, clement body loads 

>a ——- _ VECTOR, or, INITIA, read in utility vectors or initial conditions 


MESHCK > — Compute, Bandwidth, Equation blocks, block size, check Jacobian 


viscoe and implicit 


SOLVE 
assemble 
CORE IN 


a beginning sequence flag 


explicit solve 


EXPLICH 


assemble 


on tape on tape 


DISTRES 


SOLVEQ 


VELAG 


RESOLVE 
RESOLVEQ 


SOLVEQ 


RESOLVE 


—— 


DATE TTR 


Appendix E 


SAMPLE JOB CONTROL CARDS 


25 


SCHARGE eL 23) TeXrO0R. 
MONKEY oCMISU Our CL 2620000 T4000 s 1040009 THisPb. 
REQUEST FEARerIT*©® (12313) 
SK LPF (CFF ARPe lel Ten) 
Coprnk (FCAp se) Go) 

wEWwIND (FEAp) 

UNLUAN (FEAR) 

REGUF STs TAPE 30 er. (SAVED 
atTacHi1O0,QUICKIU) 
HFLO1S00006 

SET 0) 

LOAL(LGO) 

Loant1o) 

EAECUTEe 

EAT. 

exit. 
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Appendix F 


HIM CODE LISTING 


| ‘ 


SCHARGE pL 23, Terry nRe 
MONKEY 90515 He ACE AAEHOE A THUOD es LUGOCUSTH 1 PH 


HEGUEST FEAr ere (L2313) 
SKIP (RF AR ele l Tet) 
CUPYHF (FEAR. OG) 
REw IND (FE Ap) 
UNLUAD (FE OP) 
REQUEST» TARE Quer le (SAVED 
ATTACH TU, %UTCRIO) 

REL eTSO0UN. 

s€T to) 

LYALILGO) 

LOALITO) 

eReECuUTE. 

tALl. 

exLl, 


ee 


t 
PEUGRGM FEAR7T2 CINPUTS65eUUTRUTy, TAPES INPLT se TARPEABOUTPUT?® 
| © TAPE TEVETS yp TAPES, TAPE 10230755 TAPE 1 12659 TAPE 1L3=20509 
® TARE Ve=enser, TAPE 30) 
c 
CAKRK NOVE NOTE NUTE NOTE 
C Two aUDL. IMFUL FEATURES ARE INCLUDED IK IRIS VERSION 
C Le AK UPTTONAL CARD MAY KE INCLYDED IMMEDIATELY PRECEDING 
[> TRE FAiRST FEAR CaRUe THIS CARD HAS TRE ENTRY ISZ0T (15)6 
c APMIS SABLE BANGE 1S SodeLFe ISZUT eLEe®5Q00+ 
C TRE PROGKAM wILL ASSIGN THIS VALUE TO TRE SIZE OF BLANK 
Cc COMMUN ANG HEL UP TO THIS CIMTTe THIS 4ILL ALLOW TRE 
C bro kam TO Cuau AND EXECUTE IN ONLY Tre AMOUNT OF CORE REGHe 
f IF OMETTEU, A DEFAULT VaLUg UF S800 IS ASSIuNED TO 1Szol. 
C ge & } ENTERED IN CC 79 UF THE FEAaP CaRD INUICATES THaT 
Cc TRE WODE KESFEGUENCING CARDS PUNCHED BY THE FEAP PREPROCESSUR 
t ARE Tu pe INPUT IMMEDIATELY FOLLUWING TRE FEAP CARD. 
CRRA NOLE NOTE NOTE NOTE 
Cc FEA 4C 
Coece FINE Te ELEMENT ASSEMRLY PROGRAM ® STORAGE ALLOCATION FEa SC 
fe FFA 6 
COMMUN ZTIMINGS THICK ye NTR] pHACKCP yNRACK 
Common bTOd) 
DIMENSION THIN (6) 
CAIA LRINA PeGolUel doles sG/ 
DA;& SPE ARs gHFE Ap, 
CALL re TINEIN Cl e6e1KIN) 
CALL sl Kal 
TRI Crs. 
HACKCEZ.. 
tlele« 
Ne acnae 
T1L206 
cEle 
c FEA BC 
Ceeee Tost must COINCIUE wITh SIZE UF ARRAY UT FEA YC 
C FEA 19¢ 
CARAK 
( CHECK + Um UFTTONAL FIRST CARD 
CRAAX 
TSdvtsshodu 
READ (Selo) TFEAY 
Lu0 FoRmMA, (4) 
BHACKSrRACE § 
Te CTF E gh et. JFEAP) GU TU 60 
} RFA (e209) 1SZ01 
j 200 Formal (75) 
TFCTSZOT LT. S90) Is2vt=sc9 
{ TF C1S20T6GT.52 000) L920T=5000 
| So CALL ARFL (LOCF (ON CLS207))) 
C Fea 12C 
CPPP® LE ASU LAHKLE MAXIMUM DARD wIDTH FES )3C 
Cc FFA 14c 
MHANUSI..00 
COLL PELMI (UT ey S2ul em ang) FEA 16¢ 
CKAAR 


C PRINT #RUGHAM TIMINGS 
C NOTE Teal 'He TIMINGS ake NOT MEANINGFUL 
C TF STACK! JUnS wEHE PRUCESSED 
CARKAK 
Tk (ws lapeglheuy TISIRICE/NTRY 
TF (tak eGle0) [2SHACKCH/yRACK 
BRITE (heSuy) NIMLolRICh eT) oNHACK sHACKCP ele 
Svuo0 FCKMal ( #7 


SP ys eae 8- 


Oa by PROGRAM TimtGs 4/ 
46H NUWMER OF TUTAL TIME FUR EACH / 
eae TRL aNGUL AK TIME FACTOR, ZATION 7 
abr baClORT ZATIONS (cP) (CP) 47 
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Pann 
ear 
ear 
e 

st 
ia 
*0eCK Mar 
sun 
Cc 
Coen OS 
c 
¢ FI 
C 
c 
( 
c 
F C 
C 
n 
L 
ri 
( 
fi 
ri 
c 
c 
t 
C 
C 
Cc 
C 
C 
Cc 


naeanecancancacannan nancaccens 


week OF TOTAL TIME Fiomw FatH & 
HACK WAL TIME REQUCTION £# 
HEULCTIUNS (cP) (CR) s/ 
PSek 1ne3,F20e3as7 0 
¥ 
Lye 
WOUTINE sANUAL 
# LaSTkUCTIONS ahi LNPUT FORMATS FOR FEAK=7q SH eeeeeernetonves 
Nive ELEMENT aNatryStS PHOGRaM 
eeahtra IS a GENERAL (FY) TNITE (CEXLEMENT (arbMal YSIS (#)X0GKaM 
wt m FuguTSmEs To THe USER MESH INpUTJUUIBPUTs ELEMENT 45SEmnLY 
a LUTION OF EWUATIONS (LINEAR, {MPLICI' ano EXPLICIT TIME 
teehoem ts WONCINE AR) > PRESCKIRED GENERALIZED MODAL FORCES > 
parr sCrine WODAL AnU EL FwEny Uta, AnD GUurrur OF THE 
f KALI Zé ISPLACEMENTS ANU FORCFESs ELEMENT MalRICES FUR Tw 
aN, 'hMee GIMENSTURAL LINEAR ELaSTict! +, SHeELLSs PLaleS, ANy 
F Joy (lap ACE Euuftiund Protrtms AXE avAjL ol Ee AL Tew ATIVELY 
ShaS May SUPPLY TRETR OwN ELEMENT CIHKARY AY PROVIDING A 
SU.HUUTINE CaLLEU ELMTNN, wRERE NN IS a [a0 nIGIl NUMgER 
POest TF yinG ree ELEmEnr SUPHOQUTINE.s EACH ELEMENT SUtRuUTINE 
Bs AT CEAST FUUK SaASIC FUNCTIONS welCh AHE UFLINEATEU KY A 
Sai toul PakameleS, [Swell THe SUBKOUTZ Nee 
EL & Pit (ig Mae NULMONOF @ NEL eNELL ONSITE WOT ZV a NVEC oO MCT eM ely kVly 
LA st ok UCR gE STEF oy VECTs 15m) 
1S ELEMENT NUMMER, 
MA IS TRE MATEXT AL NUMHERs 
wor" r i ENS TUN 2g0k 3 
WUF id wombed be ARE ee? se Be Poun PER NOUE, 
WEL TS THE NUMBER OF EXTERN, NODES PER ELEMENT 
NELE 1S DIMENSION OF ELEMENT PROPERTY aReKay, 
STF LS Tee SLZE OF THE ELEMENT SITTER NESSe 
VSIZV¥ TS thE SIZE OF UTILITY VECTORS. 
W¥EC 15 ThE NUMREW OF CTILITY VECTURS. 
MCT IS A PRINTER Link COUNTER, 
um TS A PAHAWMETES FUR wATERTAL IDENTIFICATION. 
Vitel) JS “ATERJAL PRUPERTY mMATHITX (63 CELLS)» 
AYZ(NUIMe1) AKE NODAL CUURDINATES. 
LACNEL Ded) ARE ELEMENT PROPERTIES. NOUES® ETCe 
FInnk el) ARE NWUDAL GENEHALI7ED FORCES® 
FuRCE(NSTPe2) Is ELEMENT FurRCE vECTOR to BE 
CUMPUTER® COLUMN 2@ [S LUMFED MASS 
ESTIFENSTRoNSTFD TS ELEMENT MATRIX TO BF 
CuompuTtDe 
YVECTINST2¥,1) ARE PRESCRIHEN NODAL OR ELEMENT 
WUANTITTES. TEMPERATURES ET ce 
Utnorely LS SULUTLON VECTUR, 
1S4 £5 SWITCHING PAL AMETERS 
Swale #@ MATERLAL CrHaARACTERLZaATIONe® 
TSwa2s @® CHECK ELEMENT FOR POSTTIVE SHEA © 
[Sea30 #@ ELEMENT STIFFNESS CUmMPUTATIUNG® 
[Sweas @# ELEMENT STRESSES ant) PRINTOUTEC® 
[SweSe #® ELEMEW! (OAL COMPUTATION 
[S460 #*® NONLINEAR GENERALIZED FURCES 
UTHER [Sw May HE USEU FOR SPECTalk PURPOSES, 
USEHS Can GENERATE SURFACE LOADINGS KY FPRKUVIUING SLOAN 
SUSHOUTINES (#RERE NN TS A TWO OTGIT NUMBER BETWEEN OL any 05) 
frat SPECIFY Tre LoAU HQUITNEs THE SUBROUTINE 15 ACCESSED FY 
Me Cali 'O 


SLUNN (ND LM, NDF pNUR gNPRES,IPRESSPRy xy coh S) 


#RERE IN ADDITION To QUANTITIES DEFINED ABOVE 


0 


¢ 
y 
4 
Ka 


(elelele wel sielelelelelolelelelels MelelelM@elelelal « 


acaonana 


aon 


o 


fotetoelelelelels 


rolaeelelelels Mata lel ellel ala 


aonanoaan ec 


Ste SECTION 


JNEEGRAILIJON 


FOR ELMINNS 
WOP 1S TRE UIMENSTON OF LUA 


NPRES [S NUMBE 


* 


OF LUADEUD N 


NED SURF aCe 
OVES (MAX 4) 


TPRES(8) ARE NUDE NUMHERS OF LOADED NUOUFS.« 


ROS) ARE LOAD 
NUDES. 


VALUFS AT CU 


RKESPOND ING JTRRES 


FS(HeA8) AKE TRE COMPUTED GENERALTZEN (novel) 
ForCES For EACH UFGREE af FREEDgm Ay EACH 


TPRES NUDE 


Te} FUR Data InPeul 


TApLE IS acceSSkD 


“ 


be TAILS, 


¥ THE CaLt 


CALL INTEGLILIMeNCILoNDIMsLINT STU) 


ala 


Je) ala Tyrve f 
Fach Wal 
Tyre . 

JLiUw> 


FIILt 
KEAN 
MAIER] 


NVV AL 


PUL AM 


tlemth 
GENE KD 
HLUCK 


“ 
a NI yer IVE 


12 


STUN (GeM) INTEGRATIUN FUINTS AND WEIGHTS. 


eeroTeee M MLS 
fran UR EQUAL 


STuwe 


T 


AF SET EARL 


TO LINT. 
LINT = RETURNS WITH NUMBER 
wcll 2 Q KETURNS GAUSS FUIW 


Lim s | To 5S [S NUmmER 


eCrione 


SET LIM = 


vel = 1 RETURNS pa SPECIAL 
Pele 
T 


FUR 6 


SET LIM = 2 FORK 14 Te 
VC] & 2 HETURNS TRE ANGULAR 


Set tlw = 
SET tlw 2 
SET Llm = 


VEN TLE LCaTION CaRD 


sMeNT IS PREDEDvE 
AND LIALTS ON THE 


CARD. EXCEPT aS WUTEE Treg fh 


' 
v4 
3 


S 


FOX I BT. t 
FOK 3 ple wy 
FUR 7? ele UW 


(TS ei X si eae 


Y a CAR) wrl 
aMUuNnt + vA 
i 


a 


. Tree TUENTilY CakoS MaY alSv 


leHPap TTING THe INPUT Yala CaRpS. aS SUP 


TCITLY aAND RE LARGER 


INTEGRATION POINTS. 
TS AW) WEIGHTS IN 


OF GAUSS PUINTSy)[k- 


21) GAUSS FURMULA. 
UBIC ACCURACY 
QUINTIC ACCURACY, 
INTE GRATILUN FURMULA 
tTNEAR accURacye 
UADRATTC ACCupAaCY, 
UARTIC ACCURACY. 


ch InestTle ies THe 

Ta wHIcr [MMenlatecr 
SEGMENTS May aPPe ak 
Alm TRE USER [fh 

PLIEN TrERE ake 


OLEF ERENT pala IpENTIFlcaTlon cCaRDnS. TRE SE ane 


LUFNTITyY (RESTRICTIONS) 


QSTaHT OF EacR FROHLEM (M¥sl 


AlA)e 
ChANGE OUTPUT 


# 


AGE nEALTNGS 


Comments On OurTPyTt 
MalTeRT AL CHahactTeRi zal [une 


WOMB, CARUS 


FUL AR CONVERSLUN, (PRECCUE 


BLOCK } 


ELEMENT CUORNECTION CAaRKUS® 


GEtewatTe NOUES 


IN & LIAFAR 


GENEMATE ALL MESH Datla (BUT 
FUR A pe(OR JIUTMENSTORAL KEGION wHUSE HUUNDARY 
May KE DEFINED HY 4(H)0K B(20)COLLOCATEN POINTS 
BULUNDARY CULE PRE SCHIP I TON 
Gen EXATE OR KHLUCK) 
PRESCRIRED yoUaL on ELEMENT 
WUMAL OR POLAR AND ELEMEN) 
NODAL GENERALIZEn FORCES (P 
GENERA OH BLOCK), 

SURFACE LUADINGS (SAME AS F 
ELEMENT LUADINGS (Same aS + 
CHECK CONSISTENCY OF MESH ONLY (SAMF AS SULVE) 
VLOT MESH» RESULTS (SAME AS 
INTTTAL CONYQLETON PRESCRIPTION FOR DYNAMIC 
ANALYSTS (PRECEDE HY NODALS 
COMPLETE FORMULATION AND SUL UTLON FROM FLEMENTS 
(PRECEDE BY MATERT gNUUAL UR GENERA, AND ELEMEN 


PRECEDE all U!lHFR 


RY NODAL +GENERAs UF 


PaTr HY ANY JNCHEMNT 
4H NODAL AN FLEMENT) 
(PRECEDE WY NODAL OR 
DATA (pReECEDE BY 
RECEDE wy NUnaAlL OR 


ORCED 
ORCE)« 


SOLVE) 


GENERA OR LOCK) 


MAN BAC 
MAN &7IC 
MAhw ABC 
wan A9C 
wAn 70C 
MAN 71C 
Mae 72C 
wAn 736 
Mat rec 


MAN 76¢ 
mAy 77C 
war THC 
MA 79C 
wan wo 
MAN BI 

wAr, #eC 
MAN KS 

“MAN gat 


vA 40C 
MAN GIC 
vate 9g2C 
wAy SIC 
MOAN G4 
aA ISC 
vA 96C 
war, “lt 
YAN 6 
MAN 99C 
MANY C 
wANy a IC 
ay oec 
Matt) 9 3¢ 
val 4 
Many Sr 
wanted 
Maly 7 
“a deri 
Many oy 96 
wan) P 
“MaNr 1 1¢ 
MAINT IEC 
wanr 1 3c 
Mahy1ac 
Yan} 15¢ 
want 15¢ 
MANITIC 
many pc 
MANTISC 
MAN) 2O0¢ 
MANIZIC 
wAnre2C 


MaN123C 
Man} 24¢ 

125C 
war) 266 
MAND ATC 
Ma 1 JHC 
MANY eC 
MAND 20C 
MaNIQlc 
wAny 2C 


Mall 340 
Mahl 35¢ 
MAN] 36C 
Many 37C 


Plalialalelalta alalal olall salalalalaletalslle allallelalstelalalalelel slots olelelielelelalslelelelelelel elelalallalellells(alelallelelellelelallallelelals 


Uk PMECEUC HY MATEK[ ANU BLOCK 3} 
BESOLV USE PREVIOUS PROHLE™ DESCKIPITON wIThH NFw LOAL 
ONLY (PRECEDE RY SULVE ANU NEW LOADING CAFDS) « 


PRRLIC Uvp AmMIC SULUTION KY EXPLICIT INTFGRATION® (SAME 
AS SULVE) 

IMrcLic IMPLICIT INTEGRATION OF DYNAMIC PRORLEMS 
(PRECEDE KY Save DATA aS FUR SOLVE) 

VIL oCuF WUASTASFATLC CINFAR VISCUELASTIC [NTEGHATION 
(PRECEDE t¥ SAME DATA AS FUR SOLVE) 

rOUMLE FULMIER COMPOSITION (Samet AS SOLVE) 

ADUUe ACCUMULATE FOURTER SOLUTIUN (AFTER FOURTE)D 

Stuer WORMAAL EXIT (MUST FOLLU ALL Dalla) 


e@urnute@@ee FACH TOENTIF IER IS PUNCHED STARTING IN COL 7 (LEFT 
JUSTIFIEO), 


F¥cCeESS CakUS MAY FAIST HET@CENW CACH SECTION UF DATA, FUMEVEK, 
THe Cala ty BE UStu MUST IMMEQT ATELY FOLLOW TRE TYPE CaRkuU AND 
MUST WE Tt PROPEN UKDERe NU PAKTICULAK ORUEK OF TRE TYPE 
ChHU> TS np CESSANY EXCERPT Tal Tre FEART3 CARE MUST acwaYtS HE 
TRE FIRST CAR) IN EACH SET UF DATAs ANU KESIATCTIONS MUST bE 
CHSEwvel. 


ae) POWLEM INTTIZTION gNy CUNTRUL CakyS 
Cas) ye (A*,12a6) 


coy 7 10) %2 MUST CONTAIN WORD FEAP73 
cOL 13 TO 78 OUTPUT PAGE READER 


CAR) pe (15e]%9 346) 


ce pte s NOTA = SPATIAL DIMEWSION OF PROKLEM (1 TO 3) 
c% 7 to 12 NavMeS TO te PRINTED aS JUTPUT RE AUFRS To 
cUL 13:10 18 CUORDINATES © IF BLANK SET TO 10293 aS NEEDED. 


CoO. 19 Ig 26 


CACO ae (T5901 e6An) 


cot «6g te OS NUF = NUMBER UF UNRNUMRS PER NUDE (1 FO 6) 
CO: OF “Veh FE NAMES TO BE PRINTED aS UVIPUT HEADERS UF THE 
cOL 13 Io 18 GENERALIZED ULSPLACEMERTS aNU FORCES = TF 


ete eeeee SLAanK SET TO be2e3e4e5e6 AS NECESSARY 
Cor 37 Ty 42 


CORD ge (K15¢5F 000) 


CUL 1 to § NER = MAXIMUM RUMBER OF NULES CONNECTED 10 ANY 
ELEMENT (1) 10 20), 

Co. © fo Jo NEKTRA = Increases ELEwENT mATRKIN SIZE FROM 
NUF NEN TU NUF@NEN 6 NEXT RG 

CUL 1) FO 15 TREC © COMPUTE GENERALIZEU FURCE CHECK IF 
wOr ZERO (FOR TELME TWVARLTANT ANALYSTS ONLY) 

COL té& To 2y MAN = MAAIMUM EXPECTED BANDWIDTH, VEFAULT IS 


Set TO 1006 USED aS aN ERRUH CrKECK TO PREVENT 
FUNNING @1TH &@N ORVPOUS EKROR, 

cOU 2) to 25 [Ue = AUREER Size FOR STURAGE OF HISTORY 
EFFECTS IN TIME pePENDEN! @NaLYSISs. per AULT JS 
IhUF s& psa0Tseu 


cUL 26 Tu 30 WC) = USER INTEGER CONSTANT 
cUL 31 TO 40 CURL = USEH DEFINED CONSTANT 
cO. *) Ibo 50 Cun? = USER DEFINED CONSTANT 
CO. 31 To 40 CURB(D)» LeCOURD. MULTIPLIERs LCEFAULT = 100 
OL 61 Tu 70 UNZ(2) 5 hae MULTIPLIER, EFAULT = 160 
“UL 7) To 80 UN3(3) 5 3eCOUHD, MULTIPLIER, DEFAULT & 160 


Pel) *hEMgaK® USEK CUMMENTS UN UUTPUT. (6491246) 
SUsSEWLERT CaHUS 


c% 7 To 12 MUST CONTAIN HEMARK 


ewes * ee 


MAN, IBC 
many aye 
mat lat¢ 
mA yy ai 
MANTGeC 
MANIGSC 
Manl44c 
wan a9 
MAKIG0C 
MANIGIC 
MAINT 44C 
MANIGYC 
MANISUC 
MANTSIC 


w Any S3L 
MANIS4C 
MANISSC 
vANIS6C 
wanislc 
MANISBC 
MANY SIC 
MANTAUC 
MANTAIC 
MaNy62C 
MAN TEC 
MANY AGC 
MANY 45C 
MANTAOC 
MANIATC 
MANY ARC 
MAN) 6YC 
wan) /0C 
MaNy71C 
“ANT TCC 
MaN} 73C 
many Tac 
MAN 79C 
MAN) 76C 
wAny77C 
MAN] 7AC 
MANI TIC 
MANI HOC 
MANIKIC 
wand Rel 
MAN) RIC 
MAN) a4C 
MAK S5C 
MANY BOC 
MAN)» 
man ue 
MANTHYC 
MANIGUC 
wANISIC 
MAN] GEC 
MANI93C 
MANIG4C 
MaN}g9C 
MANIGO 
MANLIT 
MAN) GHC 
MANIGY( 
MAKR?PING 


man? nec 
MAN? OSC 
wAn2nold 
wanrnTc 
MaN?PUHC 
Maney 9C 


PNOO AOAC AANNAANAAN AAE CAREC AMAL ARAN AO AACA CAO 


(oieelel etelelle/ aleleleleleleleM@alelellsialelelalal « 


ancoe 


cOL. 13 To 78 STATEMENTS TO BE OUTRUT «© USE AS MANY KEMARK MAN?Z10C 
CARDS AS VESTRED. INSERT BEFORE ANY TYPF CARD. MAN?IIC 
manelec 
Ze?) TLlILe CrRANGE ON OUTPUT (6X91246) MAN?) 3C 
wAna 4c 
eo 8% ta 42 MUST CONTAIN TITLE MAN215SC 
CU. 13 TO 78 New TITLE DESCRIPTOR MAN?IO6C 
mMANCIT 
Qea) FAECUTION TERMINATION (6X9 Aq) MAND 16 
mAn a, 9C 
cou 7 To 10 MUST CONTAIN STOP, INSERT AFTER LAST PROPLEMe MAN220C 
MAN221¢ 
MAN222C 
Je) MaTERLAL CHAR ACTERTZATION (150191206) Mane 36 
wAnare 
ecu. 4 To <S NUMMAT = NUMBER OF DIFFERENT MATER[AL CHAKACT= MAN?25C 
ERKIZATIONS 10 FOLLOW. MAN? 26C 
cOL. 7 ty l2 MUST CONTALIW WORD MaTER] MAN? 2TC 
MAN? ?HC 
THE FOLLGwING CARUS arRt SUPPLIED FOR EaCh MaTERTAL TO HE CHARAC MAN? 29C 
Text Z&U (MUST Be EXACTLY NUMMAT SETS OF CARDS) MAN23UC 
mAna gic 
CARI) ye) FLEMENT SFLECTOR CARU (1591%sASe11A6) MAND 32C 
MAN? 33C 
COL ry. Fe °S MATERTAL NUMBER () TO NUMMAT) MAND 34C 
GOL. oF to tH ELNNW = WHERE AN [S NUwBER OF ELEMENT CLASS (01 AN? 39C 
TO 46) TO WHICKH THE CHARACTERIZATION BELONGS. MAN? 30C 
Chu te ko 77 ALPHANUMERTC INFORMATION TO RE OUTPUT. MAN? 3ZTC 
MAN? 38C 
CAkM) pede FTCe #* USER DEFINED FoR EACH ELEMENT TYPE PROVIDED. MAND39C 
MAND GNC 
@e) NUDAL CaRNS (1501%,A6) MAND41C 
MAN?GEC 
ch 3 To 65 NUMNR = NUMBER OF NOUAL PUINTS MAND43C 
Coy”. 7 to f2 MUST CONTAIN NODAL MAN244C 
MAN? GSC 
SUBSEGLENT CAKUS LAST NOUAL CakD MUST NUT RE GENERATED. MAND 46C 
(Loe) Se3F 19.0) MANPGIC 
wAn ag SC 
COL } Ig 5 NODE NUMBER MAN? 49C 
cO. 1S 1 TF t OTSPLACEMENT [S SPECIFIED MAN?SUC 
cOr. 14 1 IF 2 DISPLACEMENT IS SPECIFIED MAN?SIC 
cOu 17 1 IF 3 DISPLACEMENT IS SPECIFIED MAN252C 
COL 1k 1} TF @ OITSPLACEMENT IS SPECTFIED MANS 3C 
cUL 19 1 TF 5 DISPLACEMENT IS SPECIFIED Mahacac 
cOL en 1 IF 6 DISPLACEMENT IS SPECIFIED MANDS5C 
CU. ¢1 10 39 Lb COURDTNATE Value MAN2IG6C 
COL 31 TO 40 2 COORDINATE VALUE ® AS REQUIRED MAN2STC 
Cor 41 tu 50 3 COURVINATE vaLue mAnagkc 
MAN?PSYC 
NODAL CamuS MUST BE IN UHDEN, MISSING NOYES aReE INTERPOLATED MaNag0c 
LINEKRLY FROM TNFUT ROLES* Ik SUCCEEUVING CARUS HAVE [DENT ICAL MAN261C 


BOUNDARY CONES» THIS BOUNDARY CODE WILL BE ASSIGNED TO Tre INTER] wanenec 
VENTING NULLS. IN ALL UTNER COs—S THE BOUNDARY COVE Is stt TU ZERU manos gC 


eleRe TNale UN NOUF NUMNP OR A BLANK CAKDe MANDKGC 
MAN2ASC 
el) MUN SEQUENTIAL NOUAL GENERATOR OPTIONe (150)% sha) MAN 2 66C 
MANPAIC 
CUL rte 8 NUM4ER OF NUDAL POINTS MAND BEC 
cou. 7 TO fe MUST CONTAIN GENERA MAN? BOC | 
MAN270C | 
SUBSEQLERT CAKLS (2LSeT10939F 1060) MAN?TIC H 
wAnarec 
eo + to 5 WUDE *NUMBER MAN? 73C 
cUL #& tu 10 NUDE *NUMHER@INCREME NT wHICH wILiL He SUCESSTVELY Man?74ac 
ALLED TO NUDE*NUMBER UNTIL SUM TS GREATER THAN wane75¢ 
NODE NUMBER ON FOLLUO* ING CARL (ALGERRAIC) + MAN? 76C 
cu 18 To 20 bUUNUAKY CODE*« SAME aS INPUT FUR NODAL es MAND TIC 
Ih SUCCEEUING CARDS HAVE LUENTICAL BOUNDARY MAN? THC 


CUDES»s THIS yUUNpaRY CODE WILL BE aSSIGNED TO Mad27y¢ 
THE INTERVENING NODES. IN ALL UTHER CASES THE m4nero 
BOUNDARY Cove 15 SET To ZeRO6 MANPRIC 


33 


“s 


ev . Ne eee eee eee 


a ee ce 


C COL 24 TO 3 1 COURDINATE value . MAN? REL 
C COL 3) 1G an 2 COURDINATE vaLUE ® aS RKEQUTHED MANN? RIC 
C coy 41 71 5 3 COUR TNATE VaLUF e MANZR4C 
C vAn?PAse 
¢ @lLRMINATE al Tet mL ANAK CAND © MANIRAC 
C MANPKTC 
C 4?) hKoaRY COLE FATCHeUP OPTIUN, (6X—9A6) MAN?ABC 
wn 2a9C 
c CO he a Ve ust NTAILN SUUNDA MAN?OQOC 
C VANPGIC 
c SUBSERULENT CAKUSe (at5) MANZG2C 
¢ wAan293C 
€ n tf 7 Ne NUDE NUMBER TO Rave KEDEFINED ROUNUDARY CODEe MAN?G4C ' 
¢ cO. « to 19 Nhe GENERATOR INCKEMENT TU RE BODED ALGFBRAT}- MaANh295C 
Cc CALLY 10 Neo UNTIL Sum EXCEEDS (MAX OR MIN) THE Manado 
fi ww OF Tre FOLLUWING CARU« MANDYIC 
€ cM tt Ty 15 InC (ide CLeleteesNUE) COLE FOR SPECIFYING FOKCE “AN? 9480 
C Cur 16 ty 2 Um PT sepACtmEnt BuuylARY ConVITTONss wAn2gac 
c COL eee MANIANUC 
t lsc Cl) e640 oO» FORCE SEeErETeos : MANANIC 
( I4C(1) .Gle oo VISPLACEMENT SPECIFIED, NO wAnaged 
Cc INTERVENING GENEHATIUNeG MANIDISC 
C THC) ehTe OF DISPLACEMENT SPECIFIEDS MAIWINGC 
( GENEHATE SE fwEENW MISSING NUDES IN ALGEBRAIC wAn gd 
INCREMENTS OF NKe MANAGOC 
( MANINTC 
@ lERKMIN,TE wll & BLANK Carus @ MANQOKC 
Cc manage 
t 4°3) Putarw ok CYLTNDHICAL CUUKDINaTE CUNYERSTON TO CaRTESTAN MANALUC 
€ COLRDINATES (o%en6) MAaNaALIC 
C wAnr2€ 
C CUL 7 ¥G 42 MUST CUNTAIN POLAR (Lert JUSTIFIED) MANIZL3C 
¢ MAN I2146C 
C CAWD Je Cal SeSXo2F yen) MANQ]5C 
Cc wan rpoc 
C cOL 314 TO 15 Nae INCREMENT a0DED (ALGEBRAICALLY). Ny TO N2 MANQ17C 
( cu 30 § ‘Le FIRST NUDE TO BE CONVERTED MANI1 AC 
C Co. «6 To lu Née LAST NOLE To HE ConveRTED wan ap 9c 
C cUr 21 To 30 AUe URIGIN UF PUL AR KeCOURUINATE MANIZUC 
c cOL 31 to 4u YOs ORIGIN UF FOL AK YeCCORDINATE MANIZIC 
e wAnazed 
¢ Se) ELEMENT CARDS C15,)4%906) MANI23C 
Cc MANI24C 
¢ cOL 1 fo 6§ WUMEL = NUMBER OF ELEMENTS MANY25C 
c cH, Ff ig 2 MUST CONTAIN ELEMEN mAn326C 
c MANIAZTC 
t SUBSEQUENT CakUS (alSe2Q13/2ol4) MAN 3 2nC 
¢ MAN I29C 
c CAWD |e MOANI30C 
c MAaN@Il¢ 
c CH of Torx ELEMENT NUMBER MANIAC 
Cc cO. 4&4 Ig 10 MATEHIAL NUMHBEK MAN AI3C 
€ cUr tt +O 18 NUMBER OF SUBSEQUENT ELEMENTS USING Same MANA 
t STIFFNESS maTRIX ® SAVES KECOMPUTATION OF MANQ35 
¢ SImMTLAR MATHICES, ELEMENT MUST ALSO FAVE mANAqO® 
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IF 7ZERW UNLY First ELEwENT IS COUNTED. 


INCREMENT 
too Hepe al 


To 


LENEKATUR, 


LEF AULT = 1 


UNIT Ti gNUMSTR CARDS HAVE KEEN READ 


POPE OOOE MEOH EEE DY EMEHOUEHEGODH THEME EUEFO UHHH UEOOHUSH HOH HOOK OeUHPOHEHBOUOD 


VAN EHOC 
MONERIC 
MAN BEDC 
VU ANARYC 
MANeOOC 
MWANAYGIC 
De Nie ae Pa 
YVAKHGIC 
MAR BQ4C 
whan nod 
MARAGOC 
MANAGIC 
MALAGHC 
wARAQIC 
MONTHOC 
MANTOIC 
wan nec 
MANTO3C 
VANTOHOC 
wANT ASC 
MANNII NOL 
MANTO7 
vAnT esc 
VAKTOAYIC 
MANTIUC 
MANTTIC 
van Ty2eC 
MANK?7]3C 
MANTIOC 
hr 71 5C 
MANTLAC 
MANTIIC 
w 4,7) 8C 
MANT19C 
MANT2ZOC 
wAnTaic 
“ANT22C 
MAN723C 


SLBROUTINE FELMI (UT eI SZUT @MbAND) 


MALIN ASSFMs CY PROGRAM ANU STURAGE ALLOCATION FOR VARTAHLES 


RF AL LABL 
LOGICaL 
OTMENSTUR 


OIMENSION WEWHUF ¢) 


COMMUN 


/TARES/ 


1TrS 


CRECK FLAG eg NFL 
Conwd bea) oFLAGET) oTETLE Cp2) op TYME CT) (UT LISZDT) oRF (2) 


S) 
slTPo 


COMMUANSFURTERSNFE (SC a tNoF Lo F Cork seh G 
COMMURZLALSZ LIM®SGAUSS (5465) ewGAUSS (595) 


COMMUN /LABELSZ LAbL (6) eXHE N03) eX oF HEU (6) oF RO UHED (6) »UH ee RHEL (A) oR 


KWOK eA RD 2 pAaOhUG HEAD (12) eSTART CEASE IPGeNSTReWOR (30) 
COMMONS SHAR, XIACeSPAPF (44920) 0SG (499) 9SK (993) 0% (9920) LD (120) 
COMMUNSITILE SZ XTIL (3) eUTTL (6) 
COMMOn/VELUESZ Whe Cunt eCunesCuNn3 (3) 


COMMUNSUYNawO/TIMt op NSIG (7) NI gWSTEP SOS ¢NUMPLT pNEDATA (2003) »NPRy NPL 


CUMMONRJELUTS/ 


IPLOS 


CARARKK 

€ PuT IN 

CKAKEKE 
COMMUN 


/MINeWS 


NEw (L500) 


CuM4ON TO ROLL NODE RENUMBERING 


ARKAY 


x 


CATA CONWL/EHMATERT se 6R NODAL 
EMRUUND A sh RFORCE 


PEHPULAR 


e6hEL 


EMtNe» 6HSULVE 


e6HRESULVe 


FEL 


ic 
ec 
3c 


4c 
sc 

c 
ie 


BC 

9C 
10C 
lic 
12€ 
13¢ 
14¢ 
isc 
16C 


KR BMRUGR LE bye REMARN SOM INITI As ORGENE Mae eHITTILe 
X BHAULGP sbhyISCOt sOHIMPLICe Or NtwrON/ 

7+ BLANK 6h 
OATA WEL AGORLISTS7 626/04) oM29M39MG MS OMbOMASTHE SoTL TASG/S 


c 


pala mEsarNY! 


UATM pHars yn) 
KEP INU Jo 


94H 


/ 


/ 


C#eee Int tlatlt SeyRor FUR STAKT OF PROBLEM 
Lud READ CLipSsyou0) CCrHE ADs mBw 


40 


POMBLUANS pORELUAUS s6NVEC TOR 6NME SH 


PHNKLOCK 


eOMPLOl 4 
sHMEKPLICs 


FEL 


FEL 
FEL 
FEL 
FEL 


FEL 
FEL 


17C 


19C 
20c 
2ic 
2ec 


23c 
24c 


= 


POR n 


a 


_ 


120 


lee 
CREAR 


Gc oF 


c 


CAA 


Coeeetn 


240 
Su0 


ous 
60) 


1590 


153 


CCErE aly) 

IF CCC.FU6START) GU TU 246 

LF (CCE Ge Grr EAPT3) GO TU 240 
TF (CCoFGeCEASE) RETURN 

60 TU 1V’o 

CONTINUE 

OO lee lelei? 
HEANCLISELIQE CL) 


CC 79 OF FEAP73 CARL IS NON|@ZERO THEN KEAU NUDE 
CORRE SHURDENCE ARRAY PUNCHED IN PREPROCESSORK 


pa 6 p=! ,1500 


nEw(Ipal 

IF (Mw eNE ome) GU TU 65G 
JEL 

Kel 

KEAN) (SeA)L) Lene aHUF 
FORMAL (1619) 

JtJel 


TF (Lei be959g) GO TO 650 
TF (Lected) GO TO O16 
RTE (Meee) 

STtuP 1 


FORMAL (GQunye@® FRKOK @@ NODE KESEQUENCING CaRp OUT OF ORYERe 


Co 615 T2telS 


Keane] 
IF (NEw RUF CL) eNE eG) NEW(R) ENE WRUF CT) 
CONTIAve 


GO Tu 4¢ 


READ (TIRSe page) COXTITLE 
WolM = C¢ 

bo 1S. L = LenOtm 

CC # PETEE UES 

AREDCL) = CC 

TE CCCCE GM Conk) KREUCT) & XTIL CD) 
CONTINUE 

WEAUCATHEetoou) Cooll Ice 

nbF = CC 

= jeNUF 
Leth) 
PLANK) CC ® UFTLIT) 
cc 

cc 

ce 


aa 
T 
nue 


r 
ez 
" 


KER CETHS 6 1 gat) NtNoNhFERTRA, IMEC eMHANG HUF FNC] pCONT gCUN2,CUN3 


of Ss LelenOIM 
LF (CORnF(T) cEGeed) CUNTIT) © 1.0 
IF (Mbarel bey) Mia = Me Atal 


ISlZ = ISévt 

1OT ®& psvltye2 

ICH « 

tivEr s 

NRL = oTWLE, 

NLIteL IF 

nis | 

IeLul = 

NEITAC TSE tyme anu CyeccCr 
CALL TEVLC LE YME oud 


nltemw « 
HO LO Tap eieLaG 
FLAGCLIS.FALSFe 


STORAGE TS SET FUm @ ASTF &X WSTF ELEMENT 
Phx TMoM Side FUR GST TS 1208 THIS TS CONTRULEN BY 


ASTER gc WIR ORER © WE ATRA 


4 


wATRIK AND FORCE VECTOR 
OIMENSION 


FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
Fee 
FEL 
FEI 


FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 


FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 
FEL 


34C 
ySC 
36C 
q7C 
38C 
39C 
40C 
41c 
4eC 
43C 
44¢ 
«5C 
46C 
eae 
49C 


s0c 
sic 
Sec 
53c 
S4C 
55¢ 
56C 
67Cc 


Ste 
§9Cc 
60C 
6lc 
4e¢ 
44C 
65C 
60C 
67C 
48C 
69C 
voc 
71C 


a 


es 


220 


226 


227 


IF ME CLREC C1) »MHANeIHUF, 


“Cle le CeCOUNY 

> : *. ! 

’ Melle dele NL Gl bel ReNtENeOl.20) GO TO 729 
Stax hE CIST OF NAMES FOR A TRANSFEH ADUKESS 
bE A PSe5 el LTLE Mim 

so | i? 
= jeNtLISI 
iF oh tle QKR ) 12? 
Liste 

le . "* ant t ic 
IF é ‘ fF APT) Tu 120 

- ? . ase wre TOURS 

AMES UM a heel PUL ob LE eSUL RES es CSeFORG LO eELUSVEC Oe MES ePLUSFOU GRE My 
( Lied e€eO ee suse v2 U9e25e2609250 0255024592000 350e20001318 

345 6% phase eeLUseOUePG09EUOICHN) gu 

OE MO ELL of tes oe piuev Sel Mh ongew) 
WRITE CE IRA eeag a) CITILE CI) ) oft 2032) 
lc® 
136 1 = Bole 
me ({) = TLE CT) 
bes 

KhOT Tre “wpTeRwL AL CHaraCIPRIZallUNSe 
J 8 j 
IF (hE we shy) T 1S 
FLAG(1) 2. I RULE. 

Th (Leo if et) TOA 
NUMMA, = | 

My) = 

Me 3 * NUMMAT 


Aa OR MMAT 


é 
{ 
FimyegbeiSdntl) Gu tu Thy 
f 


t 
. 


» GO TL 100 


Inhout Tet nodal LUCATIONS ANU BOUNDARY 
IF (FLaG(7i) GU TO 2e5 
FLAG Ce) Fe lhUFe 

IF(berFen) GO TO sue 
NUMA = I 

*J = 

Me = sm 

Mo SZ pe * iy IM@NUMNE 

*O £ 6S ¢ AUF ONUMIAE 
TFimyeGletsént) Gl lu 71¢ 
IF (out eFLAGI3)) MH = MD 
Cowl twue 


IF (CC cE Me CURWUII9)) GO TO 227 
S] e 

IF (CCoFUCCUORWD(I7T)) J 2 @ 

Go fo els 

IF CoWuteFLAGI2)) GO 10 73a 

CALL CONVERT (NUMNK slkOTMe OT (IMG) ) 
GO TO tes 

CONT TNt 
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i 
(Ug NUNNP eNUMEL eNUMMAT ®SNUIMeNDDE ONENONEL 1 9 
PoE ECM!) oP TOM 9) oOT OMG) oUT (MS) pL 1 (M6) ) 

. 


COnDITION CONES. 


FEL RIC 
FEL RHO 
FE, AIC 
FEL 90C 
FEL 9I1C 


FEL 92eC 
FEL SSC 
FEL 94C 
FFL 95C 


FEL 9TC 
FEL 

FEL ae 
FELynQC 
FEL LOIC 
FELYNeC 
FELIO3C 
FFL UEC 
FELINSC 
FEL) 96C 
FELIOTC 
FFL OBC 
FELIG9C 
FELIIOC 
FELIIIC 
FEL112eC 
FELI13C 
FEL YI4SC 
FELVISC 
FELII6C 
FELIITC 
FEL1)8C 
FELITI9C 
FELI20¢C 
FELI2I1C 
FEL 22eC 
FELI¢e3C 
FELY24C 
FEL129C 
FEL1260C 
FEL) 27TC 
FELI28C 
FELY29C 
FEL) 20C 
FELIBIC 
FEL1I3eC 
FELI33C 
FEL 134C 
FEL) 3SC 
FEL136C 


FeLi aze 
FEL) 38C 
FEL) 39C 


FEL140C 
FELI4IC 


FEL) 4eC 


CALL BERGEN (NUMNP ®NUMEL oNOTMONUF SNENONEL TD 9 DT (M9) eT (MG) SUT (MS) 6 
OfumMa)) 

IF (IPGeLEe9) GO TU 705 

Tp (eNUT ORL AG(3)) MU = MU © NELI®NUMEL 

FLAG( 3) = eTRLE. 

IF img.Glelsz01T) Gu tu 719 

Go Tu 1e6 


pOUNDaRY COYE UVERWHITE 


IF (.NUTeFLAG(2)) GO 10 73; 
CALL COMES (NUMNP DT (M3) 
GO 10 Ves 


InvPUT THE ELEMENT CONNECTION akRay, 


J2g 

le (FLAG(37) GO TO 215 
FLAG( 5) 20 I RUF. 

IF (LetFe®) GO TU /0¥ 

NUMEL = 1 

MA = MU 

Mo = he © WELD ®NUMEL 

GC TU 215 


INPUT THE VECTORS OF PHESCRIBED YaLUFe 


IF (NUT sFLAG(a)) GO IO 7344 
IF (LeL—eo) GO TO ZOU 

IF (FLaGlad) GO TO 247 

FLAG(4) = efRUE. 

KEADCLTPSs 1901) NolévelPICK 

NVEC = J 

Me = we 

wmA = oH « NSI2V®] 

MO = WA 

IF (Mger®leGreISZU1) CO TU 716 

CALL VECIN (NUMNE eNUMEL ONDF ONSI Ze I PICK eNVECeDT(MA) oe OT (MS) ) 
IF CTPGelLE ed) GO TU 100 

GO TU 1e8 

CONTINUE 

IF (feLF eG) GO TO 125 


INPUT BOUNDARY LOADS 


CALL BLUALS (Te hUMih OND TM Oe NDF 801 (M3) eOT (MG) 91 (MS) ) 
IF CIPG eLte 9) GU TO 705 
GO fU 1e6 


INPUT FLEMERT LUALS 


IF (LeLEs 3) GO TO 1es 

CALL ELOBUS(T eNO LO NUF sMEIeNEL LONSTE ONST ZV ENVEC ODT (MI) 9DT (M2) 6 
DEMS) oF (me) oT IMO) 6 TOMS) OUP eOT OMG) GUT (MED) 

IF (TRG eLE> 0) GO TU 795 

GO 10 1e§ 


IKPUT Tre TXT TIAL CUNDITIONS 


IF (FLAG(A)) GO TU aal 

FLAGIO) © of RUE. 

NICD 2 I 

NSICU = NUF eNUMNP 

“G = we 

MO ® mG * NICUPNFIICU 

MC s FO 

IF (mo*2® Tout eIS2U1) oO Tu 710 
CALL YECTN (NUMNP oNUMEL ONDE ONSI COs) pNICUSUT IMC) oT (MG) ) 
IF (1PGelLEeu) GO Tu 100 

Go Tu 1es 


TeO Awl THRE ULMESTUNAL PLOTTING POKTHOLE CUNTKOLs 


43 


a 


FEL) 43C 
FEL I 44C 
FEL) 45C 
FELI46C 
FELIGIC 
FFL Y4AC 
FELI149C 
FELI50C 
FELISIC 
FEL) SeC 
FEL1S53C 
FFL 154C 
FELISSC 
FELIS6C 
FELISIC 
FELIS&8C 
FFLYS9C 
FELTF0C 
FELIAIC 
FELIAeC 
FELIE3C 
FELIAeC 
FELL 45C 
FEL 1660C 
FELIYAIC 
FEL y HBC 
FFLIAIC 
FEL? 70C 
FELI?ZIC 
FFLI72C 
FELI73C 
FEL1/4C 
FEL175C 
FELI76C 
FELT?! C 
FEL178C 
FEL 79C 
FEL1ROC 
FELIHIC 
FELIReC 
FFLIRIC 
FELI BBC 
FEL1®°C 
FEL) ROC 
FELIBIC 
FEL) ABC 
FELITAYC 
FELIGOC 
FELIGIC 
FELI92C 
FFLI93C 
FELI96C 
FELI95C 
FELI96C 
FELIOTC 
FEL I9BC 
FELI99C 


FEL200C 
FELoele 
FEL Anec 
FEL ANIC 
FEL? 04C 
FEL205C 
FEL 206° 
FEL207C 
FEL208C 
FEL209C 
FEL210C 
FEL21IC 
FEL21eC 


s*ee 


350 CONTI Wut 


IF IND IMeb Gee eANLeNENSNE S4) GU TU 707 

CrECK = oTHUE. 

eRITE CT IPS sen4u) 

LO 351i hela 

IF CU NUFeFLAG(K)) CHECK = FALSE 
7} CONTINUE 

TF cemulecrecK) GO Tu 707 

KSTeEPS = I 

NVEC = & 

NSIZV = NUMEL 

TF (NDI MeFG.3) NVEC F 

TF inDimeFie3) NSIdV = 6 

IF(FLAG(4i) GO TO 35¢ 

FLAG(4) = ef RUF.e 

MH = wi 

MO = MR © NyVECENSI ZV 

IF (My .GTe [SZDT) GU TO Tl9 


Se IPLOT = NLIMw 


TOUmM = wUeMAT 
GO TO 125 


soee START THE SOLUTION FUR FIRST LUAVING CaSE 


200 IF «FLAG(S)) Go TO 4) 


ChECKSe TRUE. 

cate Fictoc Tyme si) 

0O 123 Jele3 

IF (oNUTeFLAG(Y)) CHECKe.FALSE. 


123 CONTINUE 


c 


aan 


Cc 
C 


IF (.NUTechecK) GO TO 707 
NIM = 0 
IF (CC eNE SCONWO (22) @AKD eCCeNE eC ONMD(23)) GO 
wees SET UF VISCUELASTIC/NCNLINE AR SOLUTION CALL 
NTM = 1 
NSEW = IL 
y=0 
ISIZ = IséUT = [buUF 
IOT = S172 
TF INI CDeNE eg) GO 10 515 
WRITE (11P642038) 
nICO = } 
IF (CC eE Ye CONWD (23)) NICUE3 
NSICO = NOFeNUMNP 
MG a mo 
Mo ® yp * NICDeNSICU 
IF (MO.GTeIS}Z) GO TU 710 
Do Slu vy = wGemo 
16 OT(y) B ne 
1S “7 = #0 
w7A = mv 
IF (101 eGEemM7) M7? & JOT « 1 
mU = w7 
IO} @ IbT ¢ } 
MO ® m7? * NUF®NUMNP 
IF (M0.GT.IS1Z) GO TO 710 
MB 2 MO) 


eee SAVE MESH DATA ON TAPE FOR PLOT20/30 


IFC IPLOT.LGI.9) CALL PLUT2ZD(09KSTEPSs HEAD sNUMNP pNUMEL 9OT (M3) 9 


10 515 


K OT (MG) 8 DT (HG) DT (MB) eNVECeNSIZVeNELI] @NICOs TUUMs TOUMeNDF) 


eeee CHECK TRE MeSH FUR CONSISTENCY OF INPUT DalTA 


CALL MESHCR (NUMNP @NUMEL eNUMMAT ¢RUI Me NOF ONENONEL YO LOT OIO] oMBeMAXHAN 
K gOTeUT (m1) OT (m2) 0OT (M3) 9OT (MG) pOT (MS) 9OT (MO) »OT( MT) op 1SZANEGHY 


X TBLUKyNUEG) 
Mx = 1D! = nEGB 
ARKK 
OPEN TaPE1] GuIcKIO FILE IF OUT*OF=cCORE SOLN 


44 


FFL?) 3¢ 
FEL2)4C 
FEL215¢ 
FEL 216C 
FEL217C 
FELALBC 
FEL219C 
FEL 220C 
FEL221C 
FEL22eC 
FEL223C 
FEL 224C 
FEL225C 
FEL 22eC 
FEL227C 
FEL228C 
FEL229C 
FEL230C 
FEL231C 
FEL23¢eC 
FEL 233C 
FFL?34C 
FEL?ISC 
FEL 238C 
FELA37C 
FEL238C 
FEL 239C 
FEL2460C 
FELZe1C 
FEL 2eeC 
FEL 243C 
FEL 244C 
FEL235SC 
FEL 249C 
FEL267C 
FEL 248C 


FEL249C 
FEL250C 
resstc 
FEL 252C 


FEL254C 


PEE awe 


CARKK 


401 


IF ‘16 
1ouT 
IF (I 
CALL 
IF (cc 
IF (cc 
IF ‘cc 
IF (MA 
IF (CC 


IF (TPG 


IF (cc 


LKeGI.g) CALL WUICKIO (OLTAFE11 °991S7A) 
= | 
eNEw GO) wRITE(11P 652039) 
TICTOC( TYME se) 
eEUSCURWD (26)) NSEGST 
eEUeCOnwO (2@)) GO TO §20 
*EQ* cCONWOleg)) GO TO 415 
KRANeGTeMHAN) GO TO 740 
eEQe CONWO(12) 2ORe CC e€Ge CONWO(13)) Go To 125 
*LE* 0) GO TO 705 
eNE*CURWO(14)) GO TO 405 


¢ 
Ceeee FOURTER COMPUSITIUN 
c 


eas 
“uS 


READ ( 
welrte 
IPG ® 
NTERM 
NF (ONT 
XN 8 


FORM 
IF (FL 


LA = 
of 


TF CISL 2eseL TeMTA) 


1TP5s1002) 
(11TP6s200e) 
IPG * } 
= NTERM #1] 
ERM) @ J 


JoFloF2yF3,F4 
HEAUsIPGedoF lL oF2eFaeFa 


TRE STETFFNESS FOR THE ELEMENTS 


aG(s)) 
1 
MAXRANC1)®NDEG * NOF ®NUMNP®() 
GU 10 410 


GO TO 610 


* NTM) 


~eeee ONCORE SOLUTION IS POSSTBLE 


WRITE 
NE GH 


(11h6.2031) 
= NDEG 


T6un s 0 
Coeen CHECK TO S€e IF ILEST IS TO BE MUVED Up 


Cowes AVE FURCE 


yen 


MxamMs 
KEWIN 
WRITE 
IF im? 
OO 40 
OT(M? 
OT(m? 
“7 = 

LA 
mL 
we 
"0 
lol ® 
FLag( 
IF (NT 


CALL FORMST (NUMNP @NUMEL #NUMMAT pNUIMyNOF ONENONEL L oMB ODT e OT (ML) » 
X OT (Me) oT (v3) OT (M4) oT (M5) 8OT (M6) ¢OT OT (M9) oT (M7) 9OT IMB) © 


x OTM 
CALL 
Ire cre 


FORM 


CALL SULVEUINUMNP ®NUMEL NDF 910) Mg eMAXBANeLTAONSTF eISZASNEQBOIALKS 


KOT(LA 


CALL 
GO 10 


DYNAM 


IF (NI 
WRITE 
NICH 
NSILCO 
“G s 
Lal] = 
TF (MO 


UMNP BNO AT 

VECTOR FOR NONLINEAR SOLUTION 
oJ 

v9 

(9) (OTC LT) o LEMSomx) 
eEuem7A) GO TO 409 
ee © Fed 

weK) = OT(MTOR) 

A) = UT(MT) 

-TA 

“7 @¢J¢) 

LA * NUE GeMAABAN 
MU @ NUMNP ONDE 
MMe NOEG 


NUEG@wAXHAN «© } 
5) # eTRUE® 
MeEQel) GO Tu 520 


UleHT (LA) eNSTF OLALK oNVECONSIZV pMAXBANGITREC#NOEG) 
TIOTUC( rymEs3 
& *LEe 9) GO TY 705 


[He GLOBAL STIFFNESS AND SOLVE USING GAUSS ELIMINATION 


)e0T(MS) oOT (MU) 9D T (M7) p@OUT 9 OT (M9) oOT (MK) #@NDEG) 
LICTUC(TYME 94) 
4/2 


IC SOLUTION BY EXPLICIT INTEGRATION 


CO «GE. 3) GU TO 418 
(rTp&,2938) 
2 3 

® NOF @NUMNP 
Oo 
wo © NICD@NSICU 

eGEe [SZOT) GY TO 710 


FEL 258C 
FEL259C 
FEL 262C 
FEL263C 
FEL 264C 
FEL265C 
FEL266C 
FEL 267 
FEL268C 
FEL269C 
FEL2T7O0C 
FEL2T7IC 
FEL272C 
FEL273C 
FEL 2TeC 
FEL275C 
FEL276C 
FEL27TC 
FEL278C 
FEL279C 
FEL2EOC 
FeLoatc 
FEL2B2C 
FEL 283C 
FEL2R4C 
FEL2a5C 
FEL 286C 
FEL287C 
FEL288C 
FEL 289C 
FEL290C 
FEL2gIc 
FEL 292C 
FEL293C 
FEL 294C 
FEL 295¢C 
FEL 296C 
FEL29TC 
FEL298C 
FEL 299C 
FEL300C 
FELIOIC 
FEL30¢C 
FEL303C 
FEL30C 
FEL30°C 
FEL 306C 
FEL307TC 
FEL 308C 
FELA09C 
FEL310C 
FEL 311C 
FELAL2C 
FEL31IC 
FEL 3 ,4C 
FEL aISC 
FEL 316C 
FELAL7C 
FEL31eC 
FELAI9C 
FEL320C 
FEL Cc 
FEL see 
FEL323C 
FEL 324C 
FEL325C 
FEL326C 
FEL327C 
FEL328C 
FEL329C 


Ro 


“Sehsteseeeemnaneeemeeeenerm™ 


OO 417 JE*G_Mo FEL330C 
aiT OTt(y) & eG FELIAZIC 
*\8 MC 2 MG ¢ NSICL FEL 33eC 

4D = *C © NSIC FEL333C 

ME = M0 FEL 3340 

ME = hE ¢ NSTE FEL 339°C 

Mo = GeNSTE © MF FEL336C 

IF (My o6Te ISZ0T) GY TO 710 FELIAAC 

CALL EXPLOIT (NUMNP 6p UME @NUMMAT »NU IM NOF ONENONEL LO NSTF pNVECONSI2 V5 FEL SSBC 

XM OF EML) eb im 2) OT (M9) OT OMO) @OT EMS) eT (MG) ODT IMI) pO 12) eDT IMD) » FELA39C 

X OT CMa) eT OMG) CUT IMC) pDT (MD) OOT (ME) eOT (MF) op TYME STOUT eNSICD) FEL 3400 

Go Tro 126 FEL ae1c 
27 Jour «= I FEL342C 

CALL TICTOC( TYME 91) FELI43¢ 

FEL 3aaC 
Ceees FORM Trt NEw COAL AND OO & RESOLUTION ONLY FEL 345C 
Cc FEL 346C 

CALL HFSVEG (NUMNE oNUF eM T=] MAX ANG LSZAVNEGBOIBL Ke OT(LA) eOT (MS) 9 pet ae7¢ 

x DT Emu) 6 OT (m5) UE OMT) gOT IME) @NUF @RUMNP) FEL 446 

CALL TICTUCITYME 95) FELQ49C 

420 CONTI wut 
¢ FELASIC 
Coeve COMPULE ah) OUTPUL Tre NODAL UISPLACEMENTS AND ELEMENT STRESSES FEL AS2C 
c FELIS3C 
TF CL ReCeGlegeANUelHLA.EUeg) REWIND 7 oerdeec 
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FORMAL (Liye este apmmAs ROT BEEN INPUT UR GERERKATEU ) 


FORMAL Cle Ly (20h eae Me GHP AGe gIL4s/15—9 19M NODAL KUINTS// 


Wan al eotr toh OH eee FC 1 ke 2h6)) 
Fr 
rei May 


SLHHUGTINE ERIMAT (CHE aU oe LPG eNELMe STE CEST IF DF URCE she NT) 
Commu stake Ss [Te Sel1Pe 


DIMENWeTUr ESTTE (nb etl) oP ORCE (NT) @LOUNT) eHEAU (12) 
MHELA a (NST ROSS) eu 


Neck © (Sip «6 B49 


Nl ® 

OO 76 " = | eNSLA 

No SARL ¢ 4 

IF (Ne elste 1F) Ne # NSTF 
Mw) = j 

DO thy “MY & TeMBLA 


M2 @ wi © ¢y 

IF (*eeGleNSIF) MO @ NSTF 

wHETE CL Te ee 0) MEAD eT PGeNEL Ms (UeueNleNe) 
It = ten e ! 


bO SU | & MhemMe 


PRT 
PRT 
PRT 
prt 
PRI 
er 
pri 
PRI 
er 
PRI 
PRI 
PH 

at 
PRI 
pRI 
PRT 
PHT 
PT 
PRI 
bHI 
PRI 
PRI 
PRI 
PRI 
PRT 
PRI 
PRT 
PRT 
PRI 


ert 
PHT 
PRI 
PRT 
eke 
PRI 
PRT 
pRI 
eRT 
PRT 
PRI 
pT 
PRI 
PRI 
pT 
prt 
PRI 
pr! 
PRT 
er 
PRT 
PRT 
PRT 
PRI 
ORT 


PRI 
PRI 
pal 
PRT 
PRT 
PRI 
pRI 
PHI 
PRT 
ORT 
PRI 
prt 
Pri 
ert 
onl 
eR 


ee SR OO re 


‘30 
160 


*OECK 


Ceeee 


aol 


ace 


age 
404 


*NECK 


*#DECK 


Te = thsne 

CONT Livi 

IF (NHi ei Fe) GU TU 900 
NBT = NXT © | 


ODO AUS IT F [eA 

IPRES(I) = (PRES(I]) © INCCT) 

GO TU See 

CONTINVEt 

RE TURES 

WRITE CT PP ee09d) MOR es Slice TPRES 

Ge ro 93 

WRITE (TIP Gee od) NOK OSLO TRPRES eNKT oO INC 
[Pe * ¥ 

PE TURK 


FORMAL STAI; MENTS 
FORMAL (CLS ep eh e149 Fi oti 3 /BF 1000) 


FORMAL CIM Le [2ab ody Re eHPALE 14/75 %426HULSTRIBUTED BOUNDARY LOAUS7/ HLU THC 


Oeil « GR e@ML UAL e6 (14 —2A6)) 
Poke i (/tTeey ar f TmEwSTUNAL 945054 LOAD) 
\ 


FOKMAI (CL) us {P7EL 304) 
FORMAL (L7H @@FATAL ERROR@@ 5ASe¢1l6R LUAD CakUD ERRORS1X) 
FORMAL (opr @@F ATA ER ee NUR Sela g Be eASe AK OSHMNOMES SBISsSISK O15,» 


2th UNCRERENTS wl tee INCByehSe24m wILt CAUSE NONAL ERRGKe) 
ian 
Stvul 
SULPRULT INE SLOOL UND Lm oe NUF pNUPONPRES s TPHESe HH oAZ OFS) 
DIMENSTUN [LEHES (BH) 9h KOR) oF S (696) 
COMMUNSGALSZ LIM 9 >GuUS5 (565) #8 GAUSS (595) 
CommMunsst bps x JhCy SHALE (4620) 9556363) 9 SK (393) 0K (9920) LO (120) 
NSIDE = NPRES 7° J 
IF (NUR et Gee) NSOTUE = SPHESS4 
LIM = NSILE «+ | 
COMPUTE NURWVAL PRESSURE GENERALIZED FORCES 
Lies = CIM 
TF (NUreEGel) LIMd & I 
KUP = NDP eo 1 
bO aus IT = polity 
SS = sAUSs( TI etln) 
OS = MGauss( tel lm) 
IF (L Lm tehtoy) GO 10 a) 


SS sie 
WS = je 
comti tt 


DO 94 JJ = IolIM 

TT © sGAUSS( sso Ih) 

WT = aGaliSS(JueL IN) OS 

CALL BHI CRE (TT SS eel oe NDI Mg NPRESSASTDED 
PR = . 

OO 40¢ LF yohhRED 

PR = wth © wy (LT) @SrmArE (Gel) 

PR S&S white 

DO io L * podPRES 


DO 44% Jz enor 

Foi Joi) = Fellel) © SGERDE UI OSM IBE (401) CEN 
COMTI wit 

Nor = Nur © | 

RET, 

Fay 

StU e 

SURBUL TIPE CLP) 2 Ot Lee WUOF @NUP 9 iWPHE Se LPHESsPRoAK OFS) 
DIMENSTUN ARC Ie ae TERE SUR) @PH EK) OX (3020) OF loon) 


CATE mem deg =Ger lee Het lies mGotlet4eng de/ 

ADD MuUMdcbaraL an veEwTICay ComPUNENT STRESS LOADS 
DO HO, RBbeg 

PO 405 Lelea 

FRC Lent = FRCL ORD © PR CR oS OAC EL) OPH) Shed 
FSiPond) = FLL AerD = HRC R OL DMALL eL I OPHI 2) Shey 


a 
FAL 
SLUO3 


SLRROGT ING CLEC aba Lime Oh eNOR ei PRE Se LPHES oP AZ oF S) 
DIME MSTON LeHES (RD ohh OM) OFS (Oem) 
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BLO &2C 
RLU #&3C 
HLU #4C 
PLO 45C 
BLU AOC 
BLO 4&7C 
RLo 68C 
HLUO ADC 
RLU 79C 
pce Ti¢ 
Lu 72 
KLO 730 
RLU 7T%¢C 
BLO 75C 
BLO 76C 
ALG Tl¢ 
RLU 79C 
ALU age 
ALO RIC 
HLU Ke 
PLU 83 
RLO ROC 
KLU KSC 
SL! 1c 
SL 2c 
SLU 3C 
suv 4 
Sui 5c 
SLD . -6C 
stp 7 
Stt 8C 
Sui 96 
seb 10C 
SLb 11C 
su0 pac 
SLU 13¢ 
SLU 44C 
St 15¢ 
SLE 16C 
Stu TIC 
SLU lec 
SLU 19C 
SLU 20C 
SLY AIC 
SLU 22c 
SLL 23C 
SLD 26C 
SLU edC 
SLD Rec 
SLU e7C 
SLD 2KC 
SLL A9C 
SLO 30C 
SLY ac 
sub sec 
SLU 33C 
su ic 
SLO 2c 
SU! aC 
SUI ac 
sti 5c 
Sul 6c 
stu TC 
SLL BC 
SLI 9¢ 
Stu 1ec 
Sul 1c 
Sti ec 
ae eee 


a0 


SGAUSS (5058) ewWOAUSS (595) stu 3C 


ee 


ant emeeteesaes 


ean ALS, Ll C 
CMMON sSHAt LJAC SHALE (4420) 0 SG( 499) 0 SK (503) 0K (39029) LU L120) sth 4C 
NSIbE = NERES © ] SLi 5¢ 
LIM = WSILE « I SLD 6¢ 
C®PCO COMYUIF PLiw AL PHESSUGF GENERALIZED FORCES SLI 1 
OO 604 JJ = Pel IM SL KC 
TT? = sai SS(JseL TM) cL 9c 
®T = aGAllSScjcelL IM) SL 9c 
CALL se OR OCT Teele etd een IM eNPHESeNS TUF) stD yIC 
eR = . SLL 12C 
kA = oe Slv 13¢ 
DO She F = penPRES SLD 17°C 
RR © nH © SPAKE (4, 1Ie*% (1,7) SL 15¢ 
a2 PR = Fr © PRCT) OSmArE (ae]) SLU lec 
PR = pital ane SED: HEC 
00 406 BF = 1KPRES SLI) 18C 
BR FS SHArE (gel) OPN SLO 19C 
OO ett JEle- cLV 20C 
406 FSC sel) = FSCueld * SECA ed) @RH SE P3¢ 
404 CONT Livut SLU 22C 
ke TUR, Stn 23¢ 
EAD SLD 24¢ 

PUECKR BRAGEN 
SUPRULTINE BRKOE RN (NUMAP oe WUMELC eNO LM e NDF NEL eNELT oT COL eK YZ oF g IX) eaer 1C 
Fear LAH wen 8 862C 
DIMENSION LOUP Cp) eAYEONDIMe Do LXCREL J 97) oF CNUF 93) BK 3c 
Commun 414etsy lTe Sel TPG Rin 4C 
COMMU, “LAE CSZ LABL f6) @ XRED (3) @ XH oF HEL (6) pF He UHED (6) pUR pK 16) gRR RRR OSC 
R ghar be Aw RE Ce baGh i FenE an (le) eStaRlsCEASF es LPGeNSTRe WORLD (30) BRR AC 
COMMUN OBA, CAC SHAPE (4926) 9 SG (399) 9 SK (393) 0X (3929) LU(120) HRA OIC 
c Bun BC 
Coees GEMERLTION CF MUDgL COOKUINATES aNd BOUNDARY COVES FOR BRICK SHAPE BHR 9C 
Coeee REGIONS USTRG FLEMERF SRAPE FUNCTIONS @ ® @ BRE 10C 
¢’: Bun IC 
eeee INPUT Trt bE GION PSURERTIES akk 12C 
BRA 13C 
Bee Ci TRE epg ne) Niele eNSoh TaN eNE oA eN A Ce LHEUSE oO ITPRINT eo IMSHoLELM RRR 14C 
PEADCITHPS ei gnt) Lisidsljelayl5e16 BRR 15C 
READ CLTRS ep une) (CA C1 ed) o L193) eU=L NN) ARK 16C 
IF (NSeL Fev) NSB) Hee TIC 
TF lnkelLb ey) NE] BRA LOC 
IF ONTeLE eG) NTF) BRK 19C 
TF (wbett et) NI = 1 AkK 20C 
IF (NE Lt ou) NE F 4 KRK DIC 
IF (macLE eC) MA = J KKK 22C 
RITE CTI Re p0U9) BEAL OTR GeNRONSe NT oN ONE OMA RRR 23C 
welrE (Lieb scion) Livl2vl3e Tae lSeloe (xPEO (1) satel ey onl) Bee 24C 
OC 16 S = LeNA BRA PSC 
by WELTE CTH He ee 901) (ACT oN) oTSloNUIM) RRR 26C 
Cc Bun alc 
Cosee SET Tae cuntRUL CUNSTaNTS BRK GHC 
C BRR e9C 
NSIDE © (heel) /NDIMPe? RHR 30C 
DR @ ce e/FLUATINE) RK 1C 
OS = ce/FLUATINS) PRR 3eC 
OT = cee/FLUAT INT) KRK 33C 
NR @ yh e ) Pex 34C 
NS @ WS © ; ARK 35C 
NT santo} 4 BRA 36C 
IF ONUIMeb Gee) 4l = 1 BRK fle 
NRS = NROAS RRR QRC 
NF = WHSONE 6 NT © f HRAK 34C 
IF ONE CGE eNURNR) GU TO eng RRA 4 0C 
IF (NBC eNF ey) GO TU 30 HRA 41C 
c PRK @2C 
Coeee SET Tre HOUR ARY CONULTION COUES TO ZERO RRA 43C 
c RRR 44 
O00 20 I & Ni ohF BRK 45C 
Oo 35 5 ® LyNOF Bee 460 
15 Fluel) * org BRK 67C 
<0 Icont,) = yu ARK 4b6C 
30 CONTI nue Ree o9C 


a 
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— RL RRR ET = TURE Re — 


7}. eee oO 


eee eas 


N BNA KKK SOC 
ME 2 yh - ] BRK SIC 
T @ ft BRK S2C 
O00 304 KR * yetT RRK S3C 
§ * “fet BHR 54C 
DC 2ny J = LedS RRR SSC 
k EF = Len HRK SOC 
DO 10u IT = [ehk ARK STC 
Cay DRICK2 (ReSeTONULMeNNeNS LUE) Pan sHC 
Cc RKK 59C 
Coeee COMPUTE ThE NOpal CUURYINATES UF THE N eTR NUE BRK «40C 
. akn 6l¢ 
00 59 L = jyNDIM HER 62C 
GCs t. AkK 63C 
NO 40 M © LeNN BRR AGC 

4) = + SHAPE (Goh) #X ( Mm) RR 
Su Eyztusn) = CC = Bak Ree 
€ RRR BIC 
Coeee SET TRE SURFACE HUUNLANY CONUITIONS TU PRESCRIBED CONDITIONS BRK BAC 
€ BRK 69C 
TF (Let Qet » ICODIN) = ICOOUN) © Th BRK 70C 
IF (fet QeNk) I[COL(N) = ICOD(N) + 12 Ree TIC 
IF (jetted ) ICOD(N) = ICOD(N) ¢ 13 WRK 72C 
TF (der GeNS) ICODKN) = ICODIN) ¢ Te RRA 73C 
IF (xkernive} ) TCUDIN) = ICUDIN) ¢ [5 Rin 746C 
Tk (Ret Get T) ICUD(N) = ICOQ(N) * 16 RRR 75C 
NzNe BRAK 76C 
IE (K eFUeNT eANDe KenFel) GO To lud BRK TTC 
IF (Let @eNR .OXe JeEUe VS) GO TO 100 BRK THRO 
BE @ ME «© 4 RRR 79C 
IxX(necl oME) = MA RRR HOC 
XC pork) Fn = 1 RRK RIC 
IM( AoE) =F RRK R2C 
Tx Jor) = p © NR Ren #3C 
IXlegone) = 7 * NR © { BRA HOC 
TFINDiMeEGec) GO Iu 10 RRK ROC 
Tx (%emb) = p © WHS 1 Pun FOC 
I*X*t(hovwt) = KR © NRG * RRA RIC 
le (7e~—) 2 nH © NRO * PR HHA RHC 
Ix(FowE) = KR © NRS * RRO] Aan KIC 
log Rs R » jk akk« voc 
sie S * 3 * ps RRA OIC 
1) Tetenl REK 9eC 
IF CTRL USE 6eUe 0) GU TU 360 RRR O3C 
J 2 ([REUSE = 1)@100 6 MA RAK G40 
Do 32, 1 2 nEgMEg Reuse ARK 95C 
329 2 (NELreT) = y AHK 96C 
460 TF CLR INT oe ODLACNMEL LONE) & IXCNEL DONE) * 10 ARK OTC 
TF CIM ate F Gey) CALL PRIMSHINT ONE oNOIMe TCOUG KYL) HRA GHC 
IFC TELMebG eo 9) CALL HKTELM(NE SMES NEL gNELIO1™) BRR O9C 
TF UME .GT.NUMEL) NUMEL = ME aARRINNG 
KE TUR AkATOIC 
any WRITE (TIPE eendod BRAD OeC 
Iku = BRA TONIC 
KE TOR, HHA TO aC 
bone FORMAL (CIALS) WRAP OSC 
Loud FORMAL (OT) }) BRAT OE 
Lode FORMAL (19% g aF 1 O60) HeKDOTC 
ZOU) FORMAL CIM E  pPab a ILM eGhP AGE LOs e702 BUDE GENERATIONS // AKRY OAC 
x Pk eC SRUUMHER UF Re INCREMENTS els AKKICYC 
x PoaeeShypMsen UF S=lnCwEMENTs s1S/ Ararpoc 
A LiAeeShiC AREM UF Te INCREMENTS e157 WHAT DLE 
x LUA eCSREIRST NUDES NUMHEH e157 akn1TIe¢ 
x ToAeeSPF LOST ELFmEnT NUMBER ely/ Henry ae 
x POAeCSRELE MENT MATE 4{ AL NUMBER olS7 4X) BHA) 1 OC 
e190) FORMAL ( ann l ldo 
x PAGEL SF ACE BOUNLEHY CQUE eLios HHAT POC 
LI 1CAeeSheeF ACE BOUNDARY COVE el)os AHATLTC 
x LOA 25h feF ACE BULUNDARY COLE eT Lo7 HRA DLAC 
x TUAe2SRGeF ACE BUUNDARY COVE el yor Bway P¥C 
A LUA eeSMS5eF ACE HOUNDARY COVE elias 4AR YT 2U0C 
x VOagCSRReFACE BOUNDARY COVE eLLUs/5% 951 2a6)) arni2le 
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eUUL FORMA) (oh eye W12e5) HHA) ecC 
COD) FORMAL (Sx eG pe @P@FATAL FRROR@® [NSUFFICLENT STURAGE FUR NUDES) BRAD ?IC 
FAD HHA) 24C 
PORECK ELUADY 
SUBROUTINE ELCOADS INCU eNUTM NUP ONENONELT ONSTE oS IT LZVeNVEC eT YE sie €ta te 
Lb 1CU, AVS a PX eh eFURCE OESTIF OYECT) FLU 2C 
COMMONSLABELSS LABEL (6) oFMED CI) oAMeFHED (6) oF Me UHED (6) oURONMED (6) oR FLU 30 
Kew Awe Amn rottf OD (12) START CEASE pI PGeNSTHewORY (30) Figo 4C 
COMMUN s/SHAR, XY IAC e SHAPE (4920) 6SG (309) 0 SK (303) 0X (3929) gL OCL 20) FLU 5¢C 
COMMON 4TAKES/ [TH591TP6 a iG 
OIMENSTUN TYPE CLD eo LCOO CE) oxVSCNOLMe LoL K(NEL Lot oF INDE olde Lb 7? 
1 rURCE NOR 91) Pele) FLO «ac 
Cooee LEANINGS assuctatTen wl Tr ELEMENTS ELC 9c 
DO So. 4 = TeALN ecu loc 
MEM CL heey gag) Tt bet e INCodtEL ee FLO J1C 
INC & YawS¢yNc) €LO Lee 
TF CUNC Fey INC = } FLU 13C 
TF CJELeEweu) UFL * TEL . FLU 14C 
TF Cub eGTelg lL) GO Tu 12 FLU 1>¢ 
Ca kt FLo 16C 
SEL * VEL eLu ite 
JEL = IC FLO 1K 
teu NF Tee FLO 19C 
OC Wy JM = 1935 FLU 20C 
Tee et QewORD (UM) OO TU 2609 ELL 21¢ 
tou CONTIAWE FLO 2eC 
pe me on FLU 23¢ 
tiv WRIT) (LTH eeendh) LLU elmtyN FLO 24¢ 
Ip> =, Flu 25¢ 
KE Ture, FLO 26C 
ed NO SO, h Ss FEL eJELe INC e.0 27 
MA se we OL CLA (NEL oi) 919) Fru 2aC 
Cw & Fhe (ma) FLO 29C 
Te CoMott etbyd GO FU tto FLU 30C 
No al. 1 = yeNStTF FLO FIC 
310 FOXCE (Lo) = ye FLO 32C 
CO 33, T ® jeNEN FLO 33¢ 
KR = Lalo) ELU 34C 
Te lee, end BO TO 330 FLU 35¢ 
tu(l) = & FLO 36C 
NEL @ | FLU 37C 
OC 42: J * JeNOlh FLU 38C 
320, (dolby & AYE (eK) FLO 39C 
sdu Conttavt ELu 40C 
COLL ELML TH (Ne MAGNO IMO NUF pNEL SELL ONSITE ONSILZVeNVECOMCT oUMelteAV le FLO alc 
Y Lh ek FURLE GESTIF oP Ko vVFCTSS) FLO «ec 
Coeee ALD FuRCES py THE GENERALIZED FORCFs FLU 4a 3C 
~ REY © MET = WEL Fe FLU 44C 
TF (MCh eGTe!) GO Tu Jo” FLO 45¢ 
WRITE CLIO esn00) MEALS TPGs (FRE s Fre Je] NDF) FLO 46C 
IPG ® kG * 1 FLO 47C¢ 
wey & 5 ELU 4BC 
439 WRITE CTI bE erol) wel bePRiT) olale®) FLU 49¢ 
DO au T = poNEL FLU SUC 
Le a FLU SIC 
) IC ® (CUD(R) FLO S2c 
Ik = 190006 FLU 53C 
00 63, Jv F yeNOF FLO S4C 
. TRCIC LIL) Gu Tu 2y FLU 55C 
a4y 76 8 5G FT FLO S6C 
| TFOIC Gt elu) GO TU 464 FLU $7C 
GO 10 430 FLU 58C 
: 42u FCSeR) = FC oK) © FURCECUeT? FLU 59C 
ai Ik ® [uso FLO 60¢C 
«iv CONTE wut FLO AIC 
. #50 WHITE CT The ygone) ne fFORCE (jel) eUF1 ONDE) FLO #2C 
540 CONT wut FLO #63¢ 
6.90 CONT Lilt ELU 64C 
FETUR, FLU 45¢ 
19900) FOMMAL( I Sol te pASel ool Sesh ld, 0) FLO 66¢ 
2000 FommAy Limb e paAb BU xe eHPALE Gp 147/545) IHELEMENT LOADS//OxsemNOUEs Fro 47C 
ne ee FLU 68C 
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TT _ : 


200) FORMAL (SS ¢ 4 THELEMER To ITS e664 LUAUS OA (IT 3s] PEL 3 eGe2m 6)! FLU 4Y9C 
e€v02 FORMAT (1) ) 510661364) FLU 70C 
Zudu FORMAL (LTR wer ATAL ERROH@® 9A5925m LOAD CARD ERHUR, TyPE = 9 AG, FLO 71C 
1 lah FUR ELFMENT® [5) FLO 72eC 
ERD FLO 73C 

@OECK BCUNES 
SUMROUTINE ECQUES (NUMA, ICQU) acy \C 
aco ed 
BOUNDARY CULF GENERATOR aco 3c 
ACO 4&C 
ASSUMt IcGv ITS FILLEL ANU ONLY CRANGES AKE Tu OCCUW BCU $C 
TF NUT USER SHUULD INITIALIZE BY PLACING Tau CAROS WITH Nel ANN HCU BC 
eX =) ON FIRST CARD AND Wa NUMNP AND AX = 9 ON SECOND aco TC 
aco ec 


COMMUN/LABELS/ fOe 6 eG ree sree key ee eee 04 Od sunt Riee hei ven acl 4C 
Q 


XK AWOmD LAR GRO? AWORL A MEAD (12) oSTART CEASE s LPGeNSTHewORY (30) RCy 10C 

COMMON/SISPES/S TTP5elTP6 acu TIC 

DIMENSION [COUN (1) o ToC (6) eNHC (6) HCU lec 

h 2G BCg 13€ 

NX = | HCO 146C 

10 KP @ iy HOU 15C€ 

NYP BONA ACO 16C 

Feangiles,panuy NeN%, IRC RCO 17C 

TF (NeL EF ef ee eh oe GT eiWUMAP) GO TO S00 BCU {HC 

IT s pABs(ysC(1)) aCy 19C 

DO 20. 1 2 2e4 HCU 20C 

end TT s po*rt « [anstTeccri) RCO AIC 

ICuN(n) = If RCy 22C 

IF (NAP ee Qe 4) GO TU Loe RCU 23C 

TFC (NONE) @xPeOF ev) GO TU 300 HCU 24C 

NXP 2 = NAP HCU PSC 

WRITE CT IPA sendy) Ween RCO 26C 

300 NP = AP oe KP RCO ?7C 

TF (NKR OTe ANN eNPebt ev) GO fU Log RCO 2hC 

TF (NAr el Pe eA ON eLE ON) GO FTO LOU HCO 2YC 

iT = HCO 30C 

FREPBC OP ee teh FY 2 1 RCu JIC 

DC 404 1 = 206 acu 3eC 

LT. mm pees ACO aC 

TPerBCelyetped) It 2 It « y HCy aC 

400 CONTINUE ACU 35C 

Icon tsh) = [TF HCU %C 

GO TU Jur RCO 27C 

Cc RCU 3AaC 

Coeee OLTRPUL Tet nONZERU CONES BCU IASC 

¢ RCO 40C 

sou omcl «= BCU 41C 

Se HOU 4é&C 

BO 8). NF Le NUMNP BCU 44C 

TF CICun (hi) ep Gey) GO TO 60y HCU 46C 

veJsie} HKCU 45¢ 

neC(J) = A acy 46¢ 

| THC (J) = ECUUCN) RCO afc 

TF (yeh Tebe AnD eNeL i oNUMNP) GU TU 600 HCG, 48C 

wmcrl = wot = | BCU 49€ 

IF (MCh eGTet) GO TU 61° BCU SUC 

| WRITE CL TRO ee 0n) HEADS LRG RCy SIC 
} IPG = JFG + j RCU &eC 4 

| McT & & HCO 83 

flo WRITE CT IB ecuel) (NBC (LT) e [RCCL oLl#ley) ROU D&C 

J#9 acy Sc 

| 609 CONE LE ROU 86C 

he TUR acO S§7¢ 

[~ acu Ser 

Cosese FOKMAIS RCO Sze 

} ¢ acu AC 

1200 FORMAL (CIAL S) HBCU AIC 

CIOG FORMAL (Lib e | AAP eI Rear P ALE Gp L4/725R NUNZERU HUUNUARKY CUDFS// ACU Pec 

i J C54 eT 2RRUOR He CULE 6 32D) HKCU AIC 


2OUP FORMAL (CALL yo liteel Te gu) 
2030 FORMAL EC a9n @anON FATAL ExROK@® WODF TyCREMEWl KE TWEEN KOUFSeTSs RCY ASC 
! OM Alig lag tae UF TNCURMECT STONe SIGN CHANGEL@@) wCU ABC 


- ee 


Fro “CU &7C¢ 
®OECK VECIN 
SUBHUUPINE VEC CVUMNP SNUMEL ONDE ONSTZVeLPICK ON ePRTVeVECT) VEC lc 
Cc yEC 2c 
Cooee VECIW BLLUAS TNHUI UF JSEX YEFINED VECTORS. VEE 3t 
¢ VEC 4c 
Giees TPICK Id a CODED FKAKAMETEH AD FULLOWSs VE Ds 
Te Cn 2 oy [MeL LES VECTORS ARE ASSUCTATEU wITH RONFS ve 6 
c Teton soy, IMPLIES YectornS ake a2Suctaléu 1TH DEGREES FREEDOM VEC 7¢ 
& IetCKR = Ay IMPLItS ve cTORs AWE ASSUCIATED wITh ELEMENTS VEC as 
v 

OIMFNS TUF KR( TP ePHTV C2 eN) oe VECT (NSE ZV eN) eV aURU (4) vec 10¢ 
REAL a vec yIc 
COMMON sTareSs [Te setl Poe VEC 12¢C 
COMMUNSLABELSZ LAbL (6) eAMED (5) eo XM OF HED (6) OF Me UHED (6) PURO HMED (6) ORR VEC 156 
KAW eAWORL ee AWORU Ie HEAD (12) START CEASE DL PGeNSTHewORKU (ay) VEC 14¢ 
DATA YAURLZ.HNOUAL seRPUTNTS e GHUEGe FeyRREEUUMs6rNOe FLeomeMENTS 4 VEC 150 
DaTA KLANR ee ePRINT £4 6H + 6h WOM PRINT 7 vEC At 
C vFC ),70 
Cooee INPUT HEALER NESCHIVITONS vEC AC 
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FSTIF CL ete gte)) = ESTIF CL i*leulet) © GN®AI2? © RNPAGC FLMIVIC 
375 1 = ' © AF iat 112¢ 
38u I} = Li © Bye FLM1IL3C 
490 CONTE FLMII4C 
Coeos ( “Stn T SYMMETRY eST Lb FLMIISC 
zy { = fyANSTF FILMILOC 
K @ [«} Fumy17C 
as - io FLMT JAC 
ESTIF eked) © ESTTR CLs) © ESTIF (JeF) FEM T THC 
4s Fotth eter) s ENTIF Chey) FLY 20C 
IF (COeF Us Ze TUR FLT 21¢ 
Wy =let 
19 EctlFalel) « (lel) « CAPFURCE fC 12) 
KE “ FL) 25C 
STR&ES> MeufalyOn FUR Tre VISCORLASFIC ELEMENT FLY I126C 
vy Last) Los MESTORY COMPUTATION aNy ELEMENT ASSEMRLY FROM TP pew 7rc 
CONT £)ut FLml2sc 
MPUVATIUN AN uiPul STRESS aNnD HISTORY UF VISCUFLASTIC SOLNe FLMI29C 
Tete INTEGRATION ani IheUT The STRESS eum) 30¢ 
PG 2 ulleMad FLMIIQIC 
i ha os lewa) ELE132C 
i = UtqghAy Fury 33° 
Fe = (SeMa) FLM134C 
‘ e CKECK FUR FR STURAGF ELM) 99C 
| LL & WH (LINT#@ LI @ (NK eNG® COND EMO 1) @ND IM) 72°NG) FLY t IoC 
TELL ete LSc4) CALL VEUFFE CT TRUs LT aRetionte1S2n) FLK137C 
NL # CLINT © } FLM) 28C 
| f « LL s Vorlal FLM139C 
hohe} = Nhe FLM140C 
TF (eS, GLU) eNF et) NUMAN = el RUE FLMIOIC 
| ne = eb AL SEs Fumlecc 
4 TP (NUsRPL ToL Ee) GU IU 216 FuM143C 
' f 27> WR & J ekhUMYLT FLM144C 
; 2u5 TE CONE ALA Cie 1) cE Wee BNI) eNEVATA (NP 2) of NPL = .IRUESs €LM145¢ 
210 Conrinue Fima6l 
' LF (ibelt eid GO TO 220 FLM1e7C 
: S = U, LMI aBC 
f F @ Be FLMTa9IC 
f ! a FL™1S0C 
220 S = Siuwtlsel~)) tLMIsic 
4 1 © Slum(2ecl*1) Fimyseac 
wa © sTUwleyllol) FLMIS3IC 
230 CALL HRICR2(SeTeetoeZe WEL eNSILE) FLM154C 
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" 
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¢ 
Ir 
RNG 
A? 
3 


tr 
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UNCLASSIFIED 


Coeee 


Loo 


2u0 


140 
15u 


SS PR ee a 


140 
170 


bEA 


t 


LAr 


LOGICal PRR eNPL eNUPRNT 
DIMENSTUN GU Ce tel) oF OTTE ONSTFE ONSTE?) pFOKCE CNSTE 02) oSTUW (494) gXX(2) 6 


x Te OnFUp eb) eG ONOF OL) eUERS (302) oKIG(392) oh 1) oe SMFU(6) 
DIMENSTON VECTONSI ZV eh VEC) 
COMMON ZUYNAMOsSTIME ONSIGO7) oNEeNSTEPR SOT eNUMPLI eNEDATA(2093) gNPR NPL 
COMMON SLABREL SZ LADL (A) oe AMEDD (3) Xe FHED (6) oF Ry UHED (46) pH eRHED (6) 9 RR 
K eAwOrn be Bar Sees haOkDJerF aL (ied) oSTARTSCEASE ® 1LPGeNSTRewORD (3y) 


COMMUBsS>R AM, KIAC oe SMAPE (4420) 6SG 6353) ¥5K (393) 0X (3920) »LD(120) 
Commu sft 4eE Ss ELGS st 166 
c in 


meee ESL NUL (6029) Yi (A e20) @UUL (4920) 9 CO 


COMMUNsSE TAPES TTR LQelITPlagLIRO el I wk 
ComMUnsPELUTS/ IPL E2U 
DATA SHEDS aR RRa Br ely tee 22 93h TI ep BREAKING] pe 6RPRINC2/ 


Oat 


A 


r Te TE/ ie 2333533393333339 6060666066666 T/ 


AXISYeMETRIG FLUIU CLEMENT (UNE PUINT QUAQHATURE) 


If ({Sawestes) GO Tu 4009 
NSIHE = (NEL +1) 74 


GC 


T 


GO Clecote4eSe4)e 156 


CeaAwectTewi de THE MATERIAL FOR ILSUTROPIC LINEaR VISCOELASTIcCIITY 
BEAN CLT RPS e L090) NGoish el (49MA) 9 (99M) 9 (3 —MA) 


WREIE CT IPE edg0) Gomme (a omAd UI 5 9mA) 9U(99mA) 


CALL ANTEGL(LeQeNUIMeLINT eSTUM) 
Ch = oe 

Of = . 

OTr = ey 

Ne Ss |{ 

Oc letiad) = wE © 20) 

O(2eMA) = we * Gl 

NU @ 5 * 2enG 

It (NGeEYey) GO TO Jy¥0 
REROCATRS ei gad) (ULNA) oT e6eNU) 


BRL Debye HO1) Lit l ae Ad pe Le6 eh) 
WReE% 21) GO TO 2u0 


BEAL (Cy TeSe todd) (OL eMA) eT SNLeNU) 
wRITE CT Ib G een?) (Ol em4) sLaNncenu) 

TF (NU eL Feb 4) RE TUne 

WRITE (CL IP6 se e030) MAetiGehKeNU 

Ibo = 

HE TR 

CHECK 45h €OR NEGATIVE JaCcUnTaNS aT ELEMENT CORNERS 
CALL RIND EMenFL ol 4) 

NG = Hite) 

Ne = CcehA) 

Ne = ee (CLINT O)) eR eo (C(ND TM OT) @NUIM) S24NG) 
CON, © aS eShTANI J IS Co 

RETOUR, 


COMPUIE Tre vISCULLASTIC STIFFNESS MaTKICES FOR SHEAR aNU RULK 


SF i 
he 
NG 
GG 
RO 


Jw fee wATERTAL PROPERTIES 
4 

Chea) 

LU laeA) 

yl jeMa) 


IF(NGeLE SH) GO TO 150 


NL = * NL 

OO 14. - = eeWlLeo2d 

G6 F oG © yt loMal@mIst(ulsn (le) Mad 
KA fs (2am) 

Xn (SeMa) 

TF inrcitet) GU TO Ll 

metal *% @ 

NE FS WL & GaertK 

00 16, L 8 [lLeNlboe 

Re Ss RK © LCL OMAP PRISTIUT/DCLelema) 
Cyl = tK © FT8OG 

Ci2 = &R © [T4UG 

Cie 8 Ci? 


Cle 


ELM 
ELM 
FLM 
Fur. 


FLM 
FL 
FLM 
ELM 
Fo 
FLM 
FLM 


FLY 
FLM 


EL 
ELM 
ELK 
FLM 
-LM 
FLé 


FLM 
FLM 
ELM 
FLM 
PLN 
FLM 
FLY 
ELM 
eL™ 
er 
FLM 
FLM 
FLY 
FLM 
Fue 
FLM 
FLM 
Fur 
FLM™ 
FLY 
FLY 
FL 
FLM 
FLM 
FUR 
ELM 
FL 
FLM 
ELM 
ELM 
Fir 
FLY 
ELS 
FLM 
FLM 
FLY 
ELLY 
FL 
FLM 
Fie 
FL 
FLM 
FLM 
FL 
ELM 
Fir 
FLM 


3c 
4C 
Sc 
6c 


ae 
6c 
9 
luc 
11¢ 
12C 
13¢C 


14C 
19¢ 


16Cc 
i?c 
16C 
19¢ 
20C 
2ic 


236 
2-¢ 
25C 
26C 
etc 
24C 
29C 
30 
ca 
32C 
343C 
34C 
asc 
a6C 
37¢C 
36C 
39C 
4uc 
41c 
4eC 
43¢ 
44C 
45C 
46C 
47C 
46C 
49C 
50¢ 
S1c 
Sec 
53c 
54c 
5sc 
S6c 
S7C 
58c 
59¢ 
40C 
61c 
sec 
463C 
A4C 
45C 
46C 
467C 
68C 
agC 


C2e = CN) 
t23 * ¢i2 
€33 * Ci} 
cow Bef 


C4e = bb 
RRO = pol IN 

yTUe li gKn) 

Tos 2TUn(eyKn) 
we = 5TUwl a eK) 

CALL HRI C KA (Sel emteo scott lL eNSIUE) 

PR = «ey 
OO [4H L = jpeNFL 
BR = nh © KC pel MOrAEE Cael) 

DVUL = awn. e* JAC 
Tk. 24 
Do 38. 1T = beNEL 
COMPULF LUMBRED MASS 
TEM = ()VUL@SHAPE (a5 [) PRU 
PoRCE (I) ee) = FubCE CT) se) « TEMP 
FORCE CLL etec) = FUMLE(IT1%),2) © TEMP 
COMPUTE ELERENT SUTFRRESS 
BN = SHAPE (i gL) ZEUACPOVOL 
Zn = gthApE (ee1) ZK GAC@OVUL 
FT = SHAPE Cae T) ZKROUVOL 
Cl1@in © Clee tl « Cl4@sN 
CLAPZy © CLS Pcp, 

ClPtnr «¢ C&c#E! 

COVEN 

Claman 6 CAsett © CIgedn 

Caqefn © C34erRN 

Cla®kn © Cas#in 

C 4m eZ, © Coat, 
v= kT 
OG 475 4 = [oh€t 
BR = SHAPE ( {oud SK UAL 
2h FE SHARE (eg) JF UAL 
ET = SHAPE (ged) /RK 
ESTLFGQEL “sgt “}-= SSTIECHE “sdk 4 KN@ AL) ET*®A21 «© CN#AG) 
ESTIF GLE eytel) S ESTIF CT) ssl%l) * 2Nea31 * RNAS) 

IF (je wed) 

KESTIF CT bebo) 9 F ESTIF CT teleyl ) 
FSTIF GL Lebeouteld = ESTIF CLL elssle)) 
J} 2 gt © NUF 
Il @ {1 ¢ NUF 
CONTI WUE 
CONSTRUCT SyMMETRIE 
NO 395 1 = eeNSTF 
aK 2 [-] 

UC 395 J = 4ok 
ESTIF (Lou) = ESTIFE Chey) 
ESTIF Gye) = ESTIF (ied 
TF (Ch eE ee) RETURN 

No Yo. [elenStTe 
FSTIFCLot) = ESTIF Gist) 
Re TORK 

STRESS > COMPUTATION FUH Tre VISCOELASTIC ELEMENT 


od 
N ~ 
aa - 
neh ne 


PP 
ee 
Ne 
nou 


° 


HNPATS? ¢ ET#A22 + LNeA4S?2 
c2N@a32 * RN*®AG2 


oe 


tSTIF 


* ESTIe (sed) 


+ CO*FORCE (Lee) 


VISCOLLAST TG LOAD HISTORY COMPUTATION ANDO ELEMENT ASSEMALY FROM TR 


CONTI ve 


COMPUTAlTUN ak) OUTFLT OF STRESS BND HISTORY OF VISCUELASTIC SOLNe 


00 THe Inte Gestion 4 GuTpUuT TRE STKESS 

NG = ull amay 

NK & Ul e—hA) 

GG F ulaghA) 

mK @ (5S eMa) 

CRECR FUR 4 STORAGE 

LL * we & CLINT El @ fun tnGe CCN EM eT eNO TM) 72 °NG) 
TF (LL eGEe LSZH) CALL VHUFF CT TRUs LI wRementhe [S2n) 
Nt # CLINT 


bO $4%u 
NORRN! 
IF (Si 
NPL = 


LE = TeNL 
@ NPw 


GILLI NE ey) NUPRNT & , TRUE. 


oF OLSEs 


94 


Fly 70C 
Fim 7IC 
CCH FAC 
FLY 73¢ 
FLY 74C 
Fie 7SC 
fur 746C 
Le 77C 
FLY JAC 
Fix 79C 
FLY BOC 
FLY RIC 
Fim #eCc 
FLY ABC 
FLM R&C 
eum RSC 
FLY HOC 
FLM KIC 
Fur #BC 
FLY 49C 
FLM 90C 
tLe 9IC 
Fle YeC 
FLY 93C 
FLM Q6C 
Fur 95C 
FLM 96C 
FLM 9IC 
FuM GBC 
FLM G9C 
FLMIN0C 
Fumiatic 
FLMI1O2C 
FLY) 93C 
FLmM1nec 
FLM10SC 
FLM1N6C 
FLMINTC 
FLM1N8C 
FLMI09C 
Fumr10c 
FLMILIC 
Fuy112ec 
FuM113c 
FLM}14C 
FL 12¢ 
FumMit6c 
Fumi 7c 
FLM118C 
FLMTI9C 
FLry20c 
FuUM121C 


FLM125C 
FLM126C 
FLMI27C 
FLM128C 
FLM129C 
FLrl30c 
FLmlile 
Fim132€ 
FLM133C 
ELM136 
FLM)35 
ELM) 36C 
ELMI37C 


FLM139C 
—LM140C 
ELM141C 
ELM142C 


un 
oO 
ow 


3io 


3350 


340 


Coeee 


4ue 


490 


Coors 


IF (NUPPLTeL Ee)? GU IY elt 
DO POS «PR & YT eNUMeLT 

IF (NEU OTA INR o 1) cE Wee ANU eNEUAT ANP 92) of UeLL) NPL ® «1 HUE 
conti VUE 

TF (LeNFeld GO TO 2e0 

S$ 2 9.4 

T = gen 

GO TU 236 

S = StUeels ped) 

T = stuwl2e tol) 

wwe ofUw (a .LLr)) 


CALL aH ICR (SeTomjpeeeeNEL eNSTUE) 
CompuiFf tre PRESENT STRAIN AND STRAIN [NCHEMENT 


Tre = Ge 

Om 33 Ds hee 

C¢ = . 

oo KR 3 yoeNEL 

CC = LC « SHADE (Gen) ox ion) 
xx(y) = CC 

NG We Jt Joe 

CO" a: is 

oo = . 

OO 31. 8 = yoNEL 

CO = LO © SRARE (TL OKI MLL (Ge) © SHAPE (eK) MUL CE gm) 


DD F Lb * ShARE CT eR POULT Gym) & ShAPE Lyon) PUUL Tom) 
Bprotisd) = NOZKJaC, 
PS (uelel) = CO/K GAC 
XIG(I es) = GGAOEPS(Jelel) 
XPotsepey)? = KIG (isu) 
Dim = tH ¢ DFeES (lel) 
Tee = Ere ¢ MEPS Cielel) 
Re a ax)? 


| . 
ple . 
OO 35, % = yeNEL 


2 
cc = SHAPE (aah) 

= tT ¢ Co#lLL (leh) 
(Fy = PET «¢ CC#MUL(iek) 
ST = Goer Taney 
Str 2s Si 
= (Tuk & ETE2S/SHHKI PM (KK/S 20 = UGS3 0) 
= XIr 
Dim = pts e, © VEISKK 

S PeT#e.y = VIRSKRSIT 
O00 324" 3. = fue 
PERStiel) = DEPStlel) | Ulett 
PERSCi*tel) = DFRO(T isl) 7200 
UPD ATe. [RE SHEAR HISTORY 
IG = 4 
IF INGLES) GO TU 4ty 
DC ao, T * 4,NG 


= UC LG 1 MA) 

CO F UCB emp eMISt (ute suny 
OW = cP *0T/OD) 

OO 432 J = {92 


wee 
rove) «© CODERS (lsh) 
HIG(JeR) & cIG(Jen) © Fe 
EE = cFE*h0 
Rlot(Keles) = KIG(Neleud) © EE 
Hire) = FE 
Ne 2 He 4 
FRE = isfy) © CC#DET 
Sh. @ sf -« €s 
FE = LFepe 
Rowe) @ Fe 
Str = Slt «+ EF 
Nr sf He ] 
CONTI wut 
UrDat, tHE BULK HISIURY 


yr 
r 
z 


"PrP rye 
a pu 

2 

& 

o 


~ 

= 

> 
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a 


as 


een eal 
z 


Ce nals aii ein oe nd 
z 
~ 
.) 


FLY I GBC 
FL MI99C 
fLeP00C 
FLYPCIC 
FLeenec 
FLM203C 
FLMANeC 
FLEansc 
FLEIAC 
FLepPnTC 
Fur?o4Sc 
FL AN9C 
FLM?LOC 
FLM71ic 
FL? ICC 
Fume 1 3c 
FLYZ LEC 


410 


4aeu 
430 


435 


TE (NAR.LE eG) 60 TU 44u 

OO aes 1 = Lehn 

1G = 16 ¢ 2 

PN = GlIG*) Mad 

FE = rm (Net) «© POLGeMaAl PH IST (UIP 
xIK = KIRK @ FE 

bE = cf Ab (ent yDU) 

YIN 3 ¥JK ¢ £€€ 

HCN) = FF 

NHS AH * ] 

CONT T wk 

pe e345 I lee 

Xfo (lel) ATG(Lodd * KIRK 
RToClL* pel) = xIGlitieol) * VIAN 
El 2 ETsAKCL) 


700) ODT 


ST = St7KAKC Ld) © ALK 


Covee OUT 
IF t 
Cl 


vUL Tee STRESSES FOR THE CURRENT Timt Slee 
NUOPRNTI GO TO 640 
=z MCT « } 


IF (MC he GT eG) 


GU TU 034 


wel le (TT PG genie) 


nEaDe TIME s LPG 


641 


IpG s UhG 6 y 
mel = 16 


WRITE CT TRG op 003) Mo (ae ( lL) oTele2d) eC (KIG¢ les) eset eo 2) ley ec) oSTe 
Jehol) eJ=ly2) 9 ley 92) FT 
TFENPC) COLE PLOATA(NUMEL SNOIMyNolLy SHEL *% 9X IG yFURCE) 


OmemAs, ((CUDERS( 


LE tpPLTZp enee 1) or TO 641 
IF (LL eNEod) GO TO O44) 
YECT Cet) ATGC1 94) 
VECT (ioe)? x1Gl2ee) 
VECT (ie 3)2 ST 

VECT (ive4) = XIG(1592) 
VECHT Oy05)2 YFeS(Ze1) 
VECT (496) = UFPS(392) 
VECT hive?) ® eT 

VFCI (wend = DEPS( 391) 
Continue 

IF ({Snetieg) GO TU Son 


Cosee COMPUTE Tre LOADS FUM Tre NEAT 


DP = awtkk 

AIG(] 91) 2 x1G(P 91) PUL 
ATOC(Lis2) © KIG (304) *U0 
AIG(A9?) = xIGl3ee) POH 
X1G(201) = xIGlive) 

Stet = (STFZRR *& VIR) @ww®XJAC 
wwe = AJACHY UMD (3 —mMA) 

1) = 

00 ofy 1 ® poKEL 

gC = SHAPE (gel) — 

OHCE (Lietey) = FURCE(L +191) 
Fe = ccewn 

OC 466. J = 402 

CO es 363} 
OC 655 8 * ive 


FORMAI(SKke574 LINEAR VISCUFLASTIC MATERTAL® AATSYMMETRIC SOLID ELE 


Tlvwe sfer 


= stFecc 


- CC wEE uly (del) 


TERKMS// 10% 9 12H GTINFINITY#e 1 PEL2Z©4/10X— 
a,lpel2.,asix) 
FORMAL (LOAS (44 SHEAR MODULUS POXeGhTIMES (LUX eLP2E} 264) ) 


450 CC = CC * SKAPE (Ry LI @XIG (eK) 
460 FORCE (Tl euey) = FURCE( ILI eel) 
47g Pe * TE * Nuk 
500 CONTI WUE 
he TUR 
Coeee ELEMENT (Uni, CUMPUTATION FOR BOY TYPE FORCES® ETC. 
L Cont Iwueé 
RETURN 
#000 [eu 2 4 
RETURN 
Ceeee FORMAL STAIEMENTS 
1000 FOKMAI(215e4F 19.0) 
1004 FORMAN CHF 169) 
2000 
‘ XMENT S/T Ly 9124 SHEAR TERMS 
| x Tineierh BULK 
K Len ReiNFINTTY Sy IPELC,4/10Kelem DENSITY 
<ool 


% 


FLM?A1SC 
FLY? 16C 
Fire ric 
FLY? IBC 
FLEAL9C 
FuM?20c 
FLY? AIC 
FLY222C 
FLM223C 
FLY? AC 
FLM?A5SC 
FL? 260 
Fumaarc 
FLY AAKC 
FLM?29C 
Fim? 30C 
FL¥231C 
FLM232eC 
Lm? 33C 
PLM? AGC 
FL™?35C 
FULM236C 
FLM?ITC 


Fine SE 


ELM?GNC 


FLM?42C 
FUN ?43c 
Fur?a6C 
£LM24SC 
ELM? 46C 
Fumra7C 
FLM246C 
FLMCSOC 
Fumasyc 
FLM?SIC 
FLM2G2 

FLMaS3 

FLM?S4C 
FLM2S55C 
ELM25oC 
FLMPSTC 
FLM25A4C 
FL 259C 
FLM?BHUC 
FLMoKIC 
ELM? KCC 
FLM P43 

FLM2A4 

FLMPASC 
FLM?66C 
FLM2ATC 
FLM2KKC 
ELM?PB9C 
FLMPTUC 
FLMa71C 
FuMprec 
FLM273C 
FLEOT4C 
FLM27SC 


2002 FORMAI (IPAs GR BULK MOUULUS 9 OAS GR TIMES (10X89 1h 2E1264)) —LM?76C 
2003 FORMAL (IMU GL 7eOP2F Leo de LPSE LS e 3/UXGASe TT 91 2X9 1PGE 126 3) FLMATIC 
2004 FORMAL Cid sp 2QOvE |, IeSe2nAeGMPAGE G1 3//5Ae 1 HMELEMENT STRESSES// FLM2THC 
1 LX yg TRELEMENT e 4X 9 OHR ORD BX 9OHZ ORD os SX eGHRR=STRESS 9 3Xy FLM?T9C 
2 GALE STRESS e IX A INLL=STRESS 9 3A 9 IMT Te STRESS/IOK sg PHRR@STRAING 3X9 FLm?2aOc 
3 GARZASTHALN ae BAe OMLLOSTHAING 3A GMT TH@STHAIN) ELM2B LC 
2030 FORMAL (GH MaTERTAL el 3s 17H EXCEEDS STORAGE so 15913h SHEAR TERMS ee15e FLMPRCC 
X 264 HULK TERMS NEF Us 4, 7H wORUS. ) FLM283C 
ERD FLMPRGC 
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99 


NK DseorNPry NEOEWN— 


13 


FeEar/ 


65a 667 «Po O7Y 652 696 Fus 7s 
TA Aas a i 79 714 eal édby 699 
664 669 449 693 TIL TNA Hay 692 
409 Te 7OR 637 679 697 B98 blo 
£42 AIn ER} 652 4b) 67> 685 687 
494K 6436 672 677 676 615 654 955 
62h 9e9 609 8e4 632 634 608 623 
fel 626 ©27 645 €12 GIS 6ie 531 
S529 Set $469 SST SAK S40 546 552 
S46q S25 532 533 $34 595 600 003 
S¥7 SUH SHO 599 SYL S592 Sel S22 
43K 459 46) SIG S69 56d 556 515 
G14 194 4642 446 S63 She SO] $54 
$53 SIP 432 437 43 Sl? Sle 519 
SH} Sa? S7h 563 SHS SBO 579) SHS 
566 Sor SUNY Sv S03 504 5SU5 446 
YS 454 493 452 GHI 429 430 43) 
347 458 457 456 45) 446 424 344 
24% 44h 47H 409) 46K 646? 6462) 426 
472 473 474 415) 47K 479) Sey 476 
4)2 4)h 47) 422 412 417? 419 426 
aly 399 40 Gul 4ne2 406 405 404 
Aww 344 $92 3609 24) 342 354 353 
362 (6) 360 359 355 339 340 Je 
WAS HA 467 3ok 373 374 JIS 376 
346 779 «Jha 357 ee S28 233 396 
330 33% 322 Bed See Bet BES FES 
305 36 3O7 Fle Bit 30 309 308 
299 P9R 29S 227 249 268 267 266 
215) PT 262) 26d e264 e271) ele 2713 
247 Ahh 2A4 269 291 792 285 299 
251 752 254 245 24h 249 25y 241 
243 «62390 2380 23T «623206 62s C225 C28 
Zin 222 220 2é0 eho 22S) jet: ei 
207 20S 488 Le? 186 385 do 177 
190 yay 192 1y3 196 197 198 194 
P53 $56 9755 56 352 36 195 166 
l2a te 26 Ee? 225 120, YEO dhe 
163 1oa 157 {Ai 7H Vv 14H 71 
bl S7 55 lel 51 50 49 101 
15 4? a7 4H (149 Te a 34 
ke 8h22 52 43 64 105 106 109 
KY 92 150 ve as 37 141 67 
53 4s a7 kes hh Ws he 67 
15) le 6 JR 142 68 24 4 
oa 147 112 #%F1S 116 al 46 Ge 
a 69 oy) je 40 4) ey 31 
lo t2 13 s 9 1 la 4s 
3650 
MATER, 
eivel hand BORESSULTER TABLE 
$3587 Smh1HUU, 0027 
eum2l FURGUS HUNES*OTPLOE 
355576 4HTO”, 269002 
FLM?I mary HONpemIh NER Tale 
SISOS, 4H /6u0, caver 
FLM22 ties, td SUL Tue =SUpAKACRNOT() 
j,2?7 o>. 2000937 
FL" 22 Sm_ aHLESS SUL LUeeHKAlN 
sed2t BUSuMve «4 G00937 
fume) NECK CUVEH 
355. “BT. 200027 
NOUAL 
t) <+226 164655 
) 2.503 “1,017 
) coh 34 =-,541 
! cell? -U0.00uU 


——_—— — <n 


100 


i 


SPAce 


30469 
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ale 


S29 
3 
Ses 


515 


645 


679 
S73 
574 


575 


od} 
593 
965 
944 
347 
548 
550 
951 
$52 
954 
555 
S56 
55R 
S59 
S69 
562 
963 
564 
5604 
SoT 
S6R 
S44 
O95 
SYK 
622 
623 
624 
65 ( 
651 


052 


S14 


123 


PE ge STE Pr ee ee 


545 549 
473 
546 S50 
476 
Se7 55) 
475 
$49 553 
4764 
$5u 554 
477 
55) SSS 
478 
$53 557 
479 
554 S558 
4Ku 
$58 559 
ah] 
557 561 
46? 
5S 562 
452 
559 563 
4ue6 
561 56s 
aus 
Sec S66 
466 
$63 S67 
467 
569 573 
4ahP 
STy ST4 
469 
s7i 57s 
690 
573 577 
491 
S?y S7R 
a92 
575 575 
493 
S77 5&1 
494 
S738 Ste 
495 
ST9 53 
a4 
58, S&S 
e97 
SBe Ske 
49a 
583 587 
499 
585 5469 
500 


99) S95 
Sos 
597 601 
S06 
596 402 
Su7 
599 603 
508 


See 


SAY 


S77 


J ee se SO eee, 


601 605 
Sug 
602 4606 
S10 
603 607 
$11 
605 609 
Si2 
606 610 
$13 
607 611 
Si4 
©09 613 
$15 
Olu #14 
$16 
611 4615 
$17 
613 617 
Sila 
ele «6le 
Sig 
615 619 
Se0 
617 #e1 
Sel 
O18 #ee 
$2 
619 623 
$23 
625 62S 
Ses 
626 630 
$25 
627 4631 
$26 
629 633 
Se? 
630 434 
SecA 
631 635 
Se9 
€33 4637 
530 
634 63K 
$31 
635 639 
$32 
637 64) 
533 
4536 442 
534 
634 643 
535 
641) 645 
536 
64e K4F 
$37 
443 647 
Sik 
645 449 
$39 
446 650 
Seu 
647 65) 
56) 
653 657 
542 
£54 658 
543 
655 459 
344 


oA) 
oh? 


643 


o/s 


orn 


Obsb 


657 


beR 


636 


631 


bee 


663 


bM4 


ee a er renee oe eee 


eh ee eee 


a 


657 BAIL Bee PRE EBS HRY OGD BHA 
Ses } 
661 665 AOE DRE CHD 69) O82 59 
' 546 1 
6H) BRS Bee PMS OY2 AGT OYR B94 
547 ’ 
68 FPR wes bRS 6964 696 6099 695 
Sen } 
ORF BRT BRE be CGE 69S 100 OYA 
5469 1 
BBS 6RG HY oOhE 6ST TOL Tue 698 
550 } 
669 64) 69, O90 TOT 703 TUG 72 
551 ! 
693 697 59r 544 7405 THX T1iv 106 
552 ) 
694 69h HY4 595 706 The TIL TOT 
$53 ! 
695 699 7H, 196 7O7 TIL The TOK 
5564 ' 
697 7E1 Tee O9# 7169 713 714 TIo 
sss ? 
TOL 763 Tus Pye 713 TIS T14K Tle 
& FOx¢CE 
4 ee GAS = e258 
12 eLuT3ay 
1 IMPLIcl! 
0.001 a i F726 I 555 1 025 
2 ? O,cl sou, 314, 
STur 


126 


Appendix H 


OTHER ASSOCIATED CODE LISTINGS 


=~ ewer - 


OP reese. 


®OECK SRMSH 


20,” 


PRUGHEAM SRMSHCINBUT®UUTRUTs TAPESSINPUT, TAPE TSUUTPUT se FAPEOZOUIPUTs 


pun cn) 
COMMU, ZJACs7 NJACe TSAeNSA 
COMMON /K/ RAC 392) 
OIMENSTUN VEN(12) 
DIMENITON KHIGO1S e159 1S) NCC) o LPT (L215) 
OIMENSTUN ATCL2U5) oTTCL215) of2T (1215) sCT (res) 


COMMUN TCLWN ELT 20010) (1290) oY (1290) 92 (1290) oH (1290) AC (1290) 


DIMEN a LUN [CURRC (6) 

INTEGeH RC 

INTEGER P 

HEAL MASS 

gata RSPR SCA LNG He OARS SAO 877 
=0,9/7 385269) 89626 
WRITE l69 2004) 


FORMAI( Yh1,45460HW EL COME 10 SKULL MESH s/s) 


KEAU (Seeng) 743 
KEAN) (Se2u7) NOIMe NEL eNSADNJACs THALF se JACUN@MASSD 
eTFact 


eu? FoRmAL(h7S) 


evle 


~~ KR 


RITE (6 e212) NUIMeNEL ONSA NJAC OLRALF s JACUN®MASSD 
olF act 
FORMAL (aun WUMBER UF DIMENSIONS eeneceesecereroreeeg [5 e/ 
4 jm NUMBER OF NUDES FUR EACH ELEMENT?® ee15e/ 
UM SUBAKACHNO]) SP SP FIED eccpeoe oeslSo/ 
20H Cone ant THICKN SS rot he Leepotee Be 
Gut HALF SHULL SPECIFIED eceseees eelSe/ 
40H JACOBLAN CHECK ONLY eecvcaveess o9TSe7 
4d MASS ULSTRIBUTION PRUFERTIESeeeeee eel5e/ 
404 FACTAL GUNES 10 BE SPECIFIEDeeeseeeeoeeyl5,//) 
REAM (Seagyy) NUMMAT 
PO Su TM = LeNUmMAT 
KE (SeS)0,) DENCIM) 
WRITE (HeSug2) IMebtN (IM) 
Conti we 
FORMAL (1S) 
FORMAL (/ecgrXetlle4) 


FORMAT (¢ 9n “ATERLAL G1559H UVENSI I YeyE)1e4e2UH (BeSEC#e2 PER IN®@4) 


HECK SACOM L ANS FUR CUNPLETELY PRESCRIBEU wESr 


vo fuse yeied215 

Tetelysy 

IF (JaCUN ep ve 0) GO TO 2 

MEAL) (SeRUuu) NPT 

Forma, (15) 

OO 2 wtFpotwyt 

hE AY) (Se9l) Not li) eo KON) Gp YONI OLIN) ® TPTON) 

Cont ue 

KEAL (SeAL09) ELEM 

DO @ RISLenwELEM 

FEAL (SeR lp) Nelms FECONNINGM) oM2) 08) 

FORMA) (ey SeseP 1a) 

ICUNN(Ne} 9) = IM 

Cow tdwue 

GO to 47 

WEBI) (Secule) SCALA® SCalYe SCALZe SCALT 

FONmA, (ark ived) 

WRITE (Heeiu4) NSICE 

FOMMAI(s 264 GRORM OF BRAIN DISCRETLZAL LON 15 eek eo ShCUREDSZ) 

wHITE (Hee iy) SCALASTAL Ye SCAL CZ ySCALT 

FONMAL (a © COORDINATE SCALE FACTOHeeeeeeeeroseet gy Se te/ 
Gum ¥ COORDINATE SCALE FACTOR eccccceecccceenl Se 30s 
gon 7 CoonOinate SCA&ce FACTO. oF Sede/ 
aur TRICKUESS SCALE FACTURescecees eeoeeeeehSe3e//) 

RTUTS(NST Cho 1) OO 304 ( (NST Eel) ese (NST ZEo1) ee3) 

NPTSNIUl 

nlel 

p4a4 

nos® 

Nest 


™ ete en a of RN 


re 


ac 


HENS] Zh Oy 
NGFENSL ZL OS 
Le tna] ot o* 
LiFN> Tee te 
LEENZTCE Oe 
L3eN51 db ey 


BLILD & HASJTC ARAY Ub 


hours 
DO it [F195 
f bik og 9 
{ ti UK Seng 
NOUt = tel 
11) SelGCiede%r= NODE 
DO 16 Lee 
{ }2 JFles 
DO 12 KFS ¢ivy 
t F = NUL ee] 
k2 
KOUE 32 Lhe 
13 ABIGlLeuen)d = wove 
{ tend 3 
14 Ww 


iN 
1S Nati Ci, edeoh) SNODE 
NTUT= f 
FOL IN Tre FOES 
NO Ze 1) %4 
IMENT Geile 
IM BURRS FROM + TO. ho AS: ft 
Jes Nyel 
Le RUNS FROM 4 TO Pa AS F 
OO Al JFSenyg 
FrRowt Tuer 


NelolseJelr=e NKIG(] 9005) 
FRONT HOT Tip 


NIG (Sede lm) = NHIG(LsueN9) 


BACK IOP 

NHLGOOsGede]) 2 NHIGCIM e995) 

Hace HO 
OL NKIGESGe seth) 2 NHIGOIP 9 UeNG) 

Dg p2 KEK Oy9 
LEFT FRUNI 

NRIGseleR) = NHIGCIL 9S eK) 
HIGHT FRURKT 

NKLGloelheK) = NBIGIL NDR) 
LEFT REAR 

NalGdvaeleny 2 NKEIGCIMe seh) 
HlGM! ME aR 


NeLOONGe TR eK) = NeIO( TP eNGyh) 
2? CONTINUE 

00 25 1#heNH 

DO 25 JFis4 

Jwanlgtles 


GUES FkRUM 


JES wILTRIn 


1 


GOFS FHUM | 
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Nil 


TO 


TO 


& 


4 


4 


an 


aoOonD 


JPENGeU 


LEFT Tor 
Nel GCL eSesie NKHIG(LeueS) 


LEFT KOTTUM 
NAIGEL ede) = NeIG( 1 eJ9NQ) 


RIGHT lor 
NeIGOLesged )= NRIG(CL »IMec) 


RIGHT KUT LOM 


NeiG(LeNGeuh)= Nei GllsJPeNg) 
25 CONtIwut 


#SuUM = 
KROUNT = § 
To BEAN (Sseloun) KheVV¥ olde | PeNQUE ern 
Loud FORMAL CSAs 3eTUsOeSeiFvttued? 
WRITE (6e1 00) KXoV¥e2Ze lh eNOUE she 
TF (nuk oF boi? Gy 10 74 
KOYN' = KOM © 1 
® (NULL) FXAPSCALK 
y(ngdt)ysyre sCALy 
2tnoderez2e SCALz 
H(ng0e) sre SCALE 
IF (NUAC ente GO) MSUN = HSUM * AINODE) 
Te} (wook sie 
GO [U 7 
2} CONT NUE 
TSp S(HSUM/RUUNT) @ee5 
IF (NAC ete J) wkITE (61001) TSA 


Lou, FOMMal (ary CONSTANT SUnahACHNULD SPACE TRICKNESS [5, FS503,2%— GHI 


XNCRESS/) 


BUILD aNy wkiip TRE 4 CUNP AKAY 


~LEM2) 

vO 33 12194 

(O 33 J85ehp 

PO 33 KsS5oNy 

IM4pei 

RECEMEMECEMel 

CALL LSTORE (NFLEMy Tyg Ko hts Ge LCONN,IMAT) 
33 CONTI] UF 

OO 37 [2 Sens 

O09 34 Je lew 

PO a4 KEQeNK 

IMaley 

NELEMENELEMOI 

COLL {STORE (NELEMe Lode Ke NHIGeLCONNS IMAT 
3@ CONTINUE 

00 35 JSS enp 

CO 35 KB) 94 

Is ATER 

NELEMSMELEMel 

CALL [STORE (NFLEMeT euek ems Ge [CONNe [MAT) 
35 Conrlaue 

LO 135 JtSen# 

OO 135 KeSenA 

Iwh Tt e5 

NELEMENELE MO) 

CALL LSTORE (MELEMs Tere h eM TGelCONNeIMald 
135 Conrliul 

DO 23> JBSynR 


130 


er —————————— a 


-_ 


Sc. tHe . 
n FMelelcenenwd]Ooe] vn eo TAT) 
é 46 T 
6 JENSenl 
Ih RSEs 
[mare | cel 
FLEM= the 
Cal LSTORE CNELEMe Loto Re NeTGe LCONN GS IMaAl 
36 Cont t 
37 CONTINUE 
OO 38 {[# Ky _nle 
{ 4 jFS ee 
DG 34 K8S5eria 
IMATEN1 3-1 
NELEMSNELEmMel 
CALL LSTORE (NFLEMo TL eueK eo NHI Ge LCUNNe IMAT) 
38 Conrinut 
2u0 FORMAL (915) 
FIND SURFACE NCDF ANU NEITGHRURSe CUDE=) FUR CU 
CALL KUSTS( S8#5969 Sth eae SeasHe 95949 LOL? 
1 Sy aa) Sy 5, oy 3, oy oy 2, Sy 
CALL FOSTS (a pSeNVe SebeND ghedeN¥ 9455,N8, 1 
} Se Se NGe Se SeNile Se SeNie Se 
CALL HUSTS(SeNGeGe GehF eGo 4ING95y GeNHOSe 1 
} Se Ne Se Se WY 3 Se Ae er Se 
CALL KUSTS (a eNGeNGe HeNDeNI08 SeNHINYe S9ND8 
! Se NFe NG, Se NIWNI Te Se NIeNI2%8 Se 
CALL POSTS RG yh eSe NB gO eSey NG eGedy NGO 04, | 
1 Ne 5 Se NGe Se 39 N9n Se Ce NG 
CALL FOSTS ING eG eNFe MA gGeNIe NIV OCON%® NIGGS 
1 “G9 Se "Ge NGe SeNlle N¥e SeNIcCv NGe 
CALL FOSTS (nh se NDp4e NGsNB ge NID NI 96e NB WND 
1 N9e NFe Se W9e NF8 39 NF%9 NF9 29 NGO 
CALL FOSTSIRGeNDONIO® NF ORReNLOe ReeNDeN1 UO 
1 NGe NGe NGe NGe NGgNI1] 9 NG NGgNI Ce NSe 
DO 46 L2698 
MeL-1 
Pele] 
CALL HUSIS(4eSeL eo See Le 4959Pe Sn4yie 4e50m 
} Ss Se i) 39 Ss Le es Ss lL» le 
60 CONT Iwue 
OC @l Leary 
Ms} 
pee] 
CALL FOSTS la aNGebe ONT Pele 4eNGers GeNBole 
1 5e NIs Le 35 We ke Ze ks be le 
4] CONTI sue 
DO 42 L=AeNE 
wtL-] 
rece] 
call FOSTS(NGeGebhe Neeeelhs NGvaehs NlOs6eLs 
1 NIe Se Ce Nile Se LoeNl2e Se Loe Ni3e 


42 CONTINUE 
00 43 Lebholie 
PELe] 
MeL-1 


Cau FOSTSINIOeNGeL® NID eNB eke NLOeNGePe NHO 
1 N9e N4o Lo Nile N9e LeNleo N¥e Le NIGs 
43 CONTI Nut 

00 44 LEAK 

MeL-) 

Pelel 

CALL FOSTS(L 4959 Meaese Legree Kegese Leary 


1 Le Se 5 
& CONTINUE 


Le 5s js Loe Se 2s Le 


Seine Gi aa oon od 


RNERS® 


i’ 
53, le 
slals 

SoN13s 
rlele 

We le 
NBey lel 
W9gN] 3» 
vlels 

Se le 
wo, lel 
SeN1L3s 
14, Jol 
NI» le 
NIM eNGe 
N9gN139 


GING OMe 
N95 Le 


NGeGeMe 
Se Le 


Noelle 
N¥e Le 


, 
Se de 


omen oe 


2 FOR ENGES® 


leNS1GelPT) 
loNolGelPT) 


LeNBIGeIPT) 
als 
LeNBIGy LPT) 


LeNSITGelPT) 
ole 
TeNb1GelPT) 
rls 
leNPIGelpT) 
NBs Lelele 
leNb1I1GsIPT) 


2eN6IGelPT) 


2eNS1GelPT) 


2eN°1GelPT) 


NLOeNGeMe 


2eNSLGelPT) 


PeNSIGel py) 


ees ee 


3 


FACESe 


FOEysnsAF FUVP SHEWWN NR 


lo 
10 


sd 
we 


NED ede Kelli eSe LeoNlieGe MONL 0059 Lethe d,e 14 
s+ Le N%e By by NVe 25 Ls N9y Le 2eNBIG,IPT) 14 


(url eSe Bek eSe Gel e6e SeMe5e Col eae 15 
Se Tn 3. 25 Lew 2s Se Le Le 2eNBIGeIPT) 15 
(ALM eb ese NVM eS5e NLOeL ehe NIDePeSe NBolL ede 16 
+ N44 Le 39 N%e Le 29 N¥e Lo 1s 2eNBIGeIPT) 16 
(Le eNGy Mee NFe LebgN10y PeGeNGe LebeNBy 17 
Se Ly Senile Le SeNi2s L® SeN139 2endIGelpy) 17 
SLL ehGeN pg Pet Geh1Ue LaNkeNloe MeNGaN1Oe LoNloeNge ly 
¥% Lt» Fonte ue N7eN]2° L* N eN1L3 2eN°1GelpT) 1 
(ysl eNGe Geb ohGe B9LeNL Ue BeMONGe gel eNAe ly 
a Se LeNlle Se Lenies Se LoN)3e 2enSIGoetPT) \? 
Sp Gelehyr* GeMeNIGe NHeoLeN1G*® ANDOPONIO® NIOeLeNAs 20 
Se We Levstie N¥e LoNy ee N90 LoNJ3e 2eNSIGelPT) ra) 
a 
ae tehe arch eke GedeKPe Ge IMehe GoUeKMy 2l 
ne 3» Jo Ry 2y dey Ky ty Jo Ke 3eNBIG,IPT) ra 
CALL POSTS (Ph Tiedothe Mtr em ee Nieto hPe NiuedPehe NIGedehMs 22 
tr MF Je no wibe je RKeN12e Je Ke NI3X Je Ke JoeNH{[GelPT) ee 
$2 CONTInut 
PO S&@ [Bhan 
LO && KBA eh 
yee i?) 
Ire Io] 
Kes Kel 
RMS Rel 
Cour FOOTS (Cp geeks IMeoS eho Teen eo TP a4 eho Let ekMy «3 
1 fe Se ne le je fhe ls ge Ke le le Ke JeNKHIGelPT) 23 


Se ContInut 
DO 55 T®heny 
MO §S ReKelin 
Ips Ie} 


Ive le} 
AP= Kel 
AME Ke} 
CALL FUSIS ET NLD ere Leg lill Ke LoN1O,KP, IMeN10gKy Te N10 9 KM, 4 
} le i Ks Tehiis Re LeNleyg Reo TeS135 Ke JoNpbIGelPT) e4 
SS CONTI Aue 
st S69 Ns 
4 =f 
Ie = 1 *1 
Ihe Le} 
JP= se 
uM? 
CALL FOUSTSEL etsqge TLeubeae IModeue [eiMeoae [Peieae 25 
i le ' Se ls ’ 4 le Je ce =) Je le FeNBiGelrT) 25 
So CONFI we 
DG S? I aber 
M1 57 4» s&,,4 
le le} 
Iv= Lol 
JRE Jel 
ve tat 
j= 
CALL ¥USTS eyONlon eoihMenlOe 1PeleNnlOs LeUMONTO® IMejeNl Oe of 
1 Is Js = ls etille ls JeN}es le JeoN130 JoeNbIGeleT ray 
$7 CONTI t 
C 
C OF ERATE ¢ RDINALTES FOR UNSPECIFIED INTERTOR WODES 
bon tiki S= 
IF (NELE™ ctw me) i nea 
f As NEBenP 


CURDS UD T))@e 
wWissamMilser 
635 Confit 
IF (NMS of ) ! 4e3 
Do 695° = WIS ‘ 
1ols 
CC H3/ KElepPT 
XT (Km) S60 
Yl(n)sie 
LV(K) Hef 
O37 CTindener 
OO Keo, TES, ,NFLEM 
| ie Ae JJ=1 9A 
j neko( jhe) 
‘ MENe J. 
Jenun(™) 
; IFCLPICS) one ) 6 GC 685 
| T#Ju7i 
Lelse 
IFla@ -« 
vene] 
Ten0 qm) 
if RED se I 
ft LeJI/4 
IF f4a% oF ue JUIRER=4 
mMSheon 
Aso (4) 
; LeJut4 
* IF (Lette zich 
b meney 
Lats i) 
ETC PERE Cpe MCPD eAt Rr ORL) 
! Yee gistlegeeVcpretendeor aq 
ZT pedlcyrelepredtnred(yy 


DO 45, KRelbehPT 

Da=cT ix) 

TF CTPLOR) .GTe © 6Oke 04 eEUe UeOIGU TO $5y 

OlE* tnd 

D222 in) 

O3=2 (pn) 

MERDSAT (KIS UG 

Yonperlimysce 

Zim RCT KIS 

OLA, (x(n el) 

O2FAm> (CY (KR) -027) 

bazaars (4K) 13) 

TF (DLeeels eGheey tli lurtel 
050 CONTI aut 

IcuuNisicuunTel 

IF (LOL eFGe OOGO Ig 663 
658 CONT int 


c 
662 CONTI NvEe 
c 
Cc REAL) IN FACTAL HONES TF SPECIF LEU 
c 
IF (}RACE.NEL}) GO TU FAY 
Cc 
c BRA) jh TRE FACIAL WODAL POINTS 
RE BI (50 500G) NTF 
welrE (Se3uuy) wer 
OO wey Lei gnPrtr 
HE a(S e9)) NeHRCOND oAGUN) oY (ND OZ(N) pIPT IN) 
wel Te (Ke SL) meh C lid onlay) o¥ Cu) OZ ON) oT eT O(N? 
681 COMTI itt 
c 
Cc hEAI) Jet FacTAL CUNWECTIVITY 
' 
- c 
| KEAD(Se3Gun) RELEMF 
; WRITE (Se 3005) NELEME 
00 wee Lele NFLEMF 
HEBD (CSe smi ty NeLCUNN( Po 1G) 6 CICONN(Ng J) putp 08) 
WHITE (Ae 4UCL) Ne TCONN (NG 19) 9 CLCONN (NG J) pe] 96) 
662 CONT Twuet 
Cc 
(- UpU ATE THE ,UMNEM OF RULES AND ELEMENTS 
c 
Ne Tsrie leur tp 
NELEMENELEMONELE ME 
GEE CO eh TWVb 
c 
: Cc IF SO SPLCIUFIEGs HALVFE | KE MESH 
c 
Tk CIRALE ef ue 3) CALL HALF CNPT*NELEM® NEY) 
(9 
C  MEFERENCE COtURDINATES 10 FRANKFOKT PLANE 
¢ 


NO fey N=LenPTt 

BIND B50 S57 eK) 

Yin) 2 @¥ bin) 

LENE ONDER STS 
669 CONT Lmeur 


é CHECK Me gh FUN NEGATIVE JACURTANS 


670 DO 209 % = LONELEM 
OO Fu> TehenFe 
KR = JCONN IN, DT) 
Th Rel ehh i = .OReAeL Tet? GO VG 
IF (Ket den) GO Tu %aS 
Kaltes) = ACK) 

- FAatAel) = TER) 

| KPdet) © ¢tK) 

9uS CONT ave 
HPSITCONN INE?) 
R281 Comntheed 


ES STR ST 


134 


— 


2 2g. 2d 


Ras] WN ON 9 4) 
RG=ICUNN (Neg) 
KS=IC e (nes) 
REs] N(Neb6) 
n? iw(Ne 7) 
KK N(New) 
NSIHDE = (NEL eLISNU IMS 82 
t#} 
IF fn Aster) L2 
Lurk = 
Evur = . 
SS. = t, 
a EE? tee 
SS = «SS 
1 = le 
oi J= {ed 
rh 2 ari 
UL = . 
O sun * bee 
Us 2 
CAL HR IAC (UUel Pe Sse NUL Me NEL eNSTDE 9X aC) 
EyuULeevur ex ac 
Cc 
C CUMPLUTRE 6ASS PROPERTIES OF SKULL 
Cc 
LUR = LURK ¢ j 
LF (x gaCeGleitev)? GUY ty 3 


WRITE (Se 204) ROK GACH CL COnNONOLL? PLL=L oWwEL? 

25K FORMAL (sER ELEMENT gle ei ky BAVACUBTANSEL1L 64/0 ICONN SRI 14) 
JAC SN 
WELIE (A eHo A) Kim Ld ACK) 0 KIKSD © KIKS) 9 KURDS) e KIKAD ® KIRK TI OK EKA) 
WHITE (Agus 4) YORL © YUR2) 9 YOR3)9 Y(RG) 9 YIKS)g YUR6)» YURT) oY ORR) 
WRITE Chew w) ZOKRL) © LiKe) eo LIKI) 9 ZIKRG)o ZUIKO) 9 Zim) eo ZIRT OL(KH) 


Ana FORMA, (ee ¥ vSELa.a) 
Bud FORMAL ler Y 6&1 aeg) 
BUG FOHMA, (Ar Z sok lo.4) 


LAR UM=EICUMRL (LUM) 
wRI ThE (fey TCONN Civ9UK UM) 
296 FORMAL (/oS~ NODE eI DSe ten 1S INVERTED) 
30 CONTI woe 
VOL = voOL «¢ EVOL 
IM = |[CONN(ne 10) 
Fréss =EyuLelEN (Im) 
MASS = wASSeF MASS 
Tk (Ma5S eb ed) CALL MASSOTS ORL eRe oR enSeKS ene Ko KHeP MASS give 
ANELEM,MASS, {HALF ) 
2u9 CONTINUE 
WRITE (6099) VULemMsS5 
G9 FoumAy (s9r VOLUME EvELLe@el OM INCRE S@@3/7hH MASS =eE Ll eae 
KR SH LHTSEC##2 PER icH) 
IF (JuCON of Ue 1) GU TO 94 
Tk (JAC oL Ten) 60 10 44 
we iE (Ae 9a?) 
942 FORMAL (755% MESH HAS HEEN CHECKEO=NO NEGATIVE JACORTANS we RE FUUN 
x0/) 
91. FoRmAy (pSejome Tie ko JF 1003915) 
00 705 NFLeKPT 
FUNCH gle Bugg th) eXEnb ev (Sh ve Ue iwe eT) 
Too WRITE (he% 1) ReBC Nd oe (ri) eV ON) OZON) oIPT OND 
PUNCH G5e (ie TCONN Give} 0) 6 (LCONN(NeM) oM2] 98) 9NE] ONELE MM) 
WRITE (eeS 4) (Ne LCONNG Ye 1G) » CLCONN (NGM) oMal ob) eNal eo NELEM) 
Fy FowmMAL (215/814) 
GC fQ 94 
940 wRITE (69441) # 
9a) FowmAy (Sk KOLEgTSecrhel I MEXCEEUS KPT) 
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CALL INCOME (DRC P eNSTH) Abs COPESTIF eUDL eMMeNSIF 92) 


r= 1 

OO «O83 1 = yeNDIM 

DO 4C3 J ® 4eNDF 

UL (JelowFL) = CCR) 
KRaRK ey 

NSTe = NEL «6 ADIM 

NL = (NUTM®(NDIMel) 72 
xn = D(4emes)y 

CO F Pe JU C2—MA) 

OO ols IT = yelINt 


CALL aRICRZ(STUmM( Le TI) pSTUW (So JIT) pSTUW( 3911) eNDIMeNEL SNSIDE)D 
CALL LnstAe stu (bell) oSOxeSMApE Cy nF od) Yim 


BO eli K = yet 

UPLIRs TI) = STG(R) 

UPL (BelleLInt) = CeStn) 

Ik (LS4ebOeg) GO TU 613 

FORM UTSPLACEMENT GRADIENTS 


OO 40> L = ysNOF 


Le 
Oc *u4% w = yeNSTP 
Teme = TEMP ¢ SHARE (Gok) MUL C1 yA) 
QOVtLes) = FpMPs x Sac 


FowM sTHESS = STRAINS 
DO ayer IT = y,NOIM 
AMT) & pe 

O00 400 J = j yNEL 


MXCL) = FACE) © ACL 9d) OSHAPE (aq) 


OLLtfelTid = xetY) 

EeS(l) = Ly(lel) 

FPS¢e) = yes (OV 92) 40V(2_1)) 
FRS(3) = Lv(ese) 

TEMy = FeS(y) + Ers(3) 

TF (NOL MeFGe ge) GO FO 408 

FPS(4) = U.S (KV (e931) 4V 6 992)) 
FRS(St = Oy (393) 

ERS(6) = yeg@ (Ovipeg) tv (301)) 
Trump = ffMp «6 €PS(S) 

bO v9 LT = pone 

Ste(l) = GueePS(1) 

00 4t 1 = fetlee 

Slbel) = SIGCT) © AneTeMp 

PU S STU lagry? 

Thee Ss wilt) (3yMA) @*R UAC 

It 4 1 

OC ale t * yoNEL 

Sk 2 share (;eT) @wu 

SY = SHAE (pel) wy 

Sh = SHAPE (ag) @TEMP 

FOKCE QI Le 1) = FukCcEeCI) ol) 
FowCE(Liebey) & FuRCE(I1¢101) 
TF (NOL 4eb Gest GU HO 4t2 

SZ = SHerE Cagl) ewu 

FoRCE tT? «¢)) FORCE II) 92) © 
FOMCE CLE * Loy) SFURCECLI%101) © 


= SIG(LI@Sx = SIG 21 eS TeSMOUNUL (Tel? 


- SIG(2) Sx 


sZ*siGle) 
SfesiGla) 


SIG 09) SV eSMPUDDL (2e1) 


FORCE (These ip sFORCE CL leeed) = SZ®S1G(5) =Sy#516(4) -SeeStG (5) 


c =smeulUL (al) 
ab S27 © NYP 

Comp taut 

NSIe 2 Anim 6 NyIM © 1 


wo *¢1 5 = tenste 
NOPRNI * NPR 


IF INS( Gt) .GTeQ) NUPKNI = STKUES 


NhL = oh aL She 
IF (NUeriT eC es) Gu lu 44 
OO af NP oe poNUmpe! 


TEOMA Oh op) cE Ue Noat ie MED ALAC NPS?) ob ed) 


is4 


NPL & TRUE 


FLM 73C 
ECM 74C 
Fie 75C 
FLY 76C 
FLY 77C 
Fume 78C 
FLY 79C 
FLE BOC 
ELM AIC 
FLY HEC 
FLe a3c 
FLP 44C 
Fie ROC 
FLM ROC 
FLY ATC 
FLY 2B 

Fim #9 

FLY OC 
ELM GIC 
Fie 92c 
FLY 93C 
FLY 94C 
Fim 95C 
FLY 96C 
FLY 97C 
ELm 98C 
FLY GYC 
FLN} OU 

Fimjol 

Fim1aec 
ELM1IC3C 
FLMJQ%C 
ELM) OSC 
FLNI9EC 
Furrolc 
FL] OBC 
FLMIC9C 
Fury ,0C 
FLM) TIC 
ELMI1eC 
Fim, 13C 
FLM)14C 
FLMIISC 
Fumi1ec 
FLEII7IC 
FLMITI8C 
Fury, 9C 
FLM)2NC 
piMmidle 
Fimye2ec 
FLM123C 
ELM) 20C 
Fim) 25C 
FLM)26C 
ebhre7c 
FLM)28C 
FLM)29C 
FLMI IC 
FLKy IC 
LE) 3I2C 
FLM133C 
FLH) yee 
tLy135C 
FLMT IOC 
eLM137¢ 
Fem) 38¢ 
Fumr a7 
ELE] 40C 
FLMYOIC 
FLM)4eC 
ELM1463C 
FLE)aec 


Crtee 


418 


419 


420 
42) 
5 


Ceces 


Tb liek TedanieeNUlelthL) GO TU wel 


OO #i4 [| S&S [eNDEM 

ERECT) & wes 

CO als 1 * yodl 

STGtl) . 

EPStl) 2 de¢ 

DO aif Lb petint 

Et (pe*Sh (Gd FSTUm Ed ehd IPC De *S eld sS Uwleell i sael 
Trin OimeFueg) © w E@ep cond lads gtuml aot) 720 


ai i 19n0 iM 


RRC 1) = RACT) © EH ll eol) 
O 427 L *® yent 
SIel]) = Ss (y) * E*uCL be) 
EPSCL) = ERStT) © E*®UNLITeLeLINT) 
i t 


O eeu 
COMPUTE ThvpRTANTS FUR T¥0 (2 PROWL LEMS 

FINDiMeME ec) GU 10 418 
Bip SIRES (ND IMesTOeS10(4) 9SIG(3) 651G1(6)) 

CALL KSTRES (NDI Meth ooF PS (a) st PS(5) ERS (6) 

“CT & MGT = 4 

IF (eCiebTesy GO Tu 419 

WELTE CT hey c0u2) Us ADs TEMES TPGe (XHED(1) @AMoleloNIM) » 
KX (Tet Cl) See TS lene) 

WEITTE CT EPG ee gal) COLANK oMLANR OLE] oOIM) o CrMEU(T) @EMe l= oNL) 
IF UND EM ete 2) WRITE CT TRA, 2006) 

Iky = JG » } 

wc} = 

WRITE CE LES (ARI) No (XE CLI o TST eNO Ih) © SIG 

WHITE (LIRA, (ARC) UMemagEPS 

TR OMRL) CALL PLOATAIN SNDIMsNe ue THEN XA eI G oF URCE) 
Cont liut 

Re TURE 


FORMAL STale ments 


1000 FONmA) (4610 ,90110) 


2000 
ev0l 
e90e¢ 
2003 


204 


Fokmal (Ske5ah LINEAR ISUTROPIC ELASTIC [mO/THREF DIMENSTONAL ELEMEN 
LTA 7 54g AME SES eS eS eGHNU BoP Lye d~SKySHEMO Heh 1D eS/S%gGHAA Sof 105 
2e2XSdr (hed FOR PLaNé STHESS, 2,0 FUR PLaNe STRAIN AND 30)75%) 
FoRmAy (x04 ¢248)) 

FORMAL (Ayo |} eAB oF D3 eS 01 TA oe OHRPAGE SIL 3S/5Ag 1 GRELEMENT STRESSES y/ 

1 L*e7rFLEM eh Te9(246)) 

FORMA) (Sx e5epHOPFALAL FRHUK 1 2ee@ ATTEMPT TO USE INCOMPATISLE MODE t 
PLEMENT wttts WRONG NUMBER OF CONNECTED NODES) 

FORMAL CTAR s pH LS VALUES 64% pHHET VALUES) IX eo SHANGLE 5 

Fr 

SUBROUTINE VSTIF (ip IM NF pNELYNSI DE yST, FORCE ySTUMW,LINT) 


COPPerevee VSTIF BPEWHTITEN@ ORO 7T/2Us [JOO COPORHOOEH OD EOEO HOH DHEMOUOO HOO HHS 


50 


St) ~ eweee 8  eee ewe 


100 
150 
200 


DIME WSTON STUW(49e7) oSTC1) FORCE (NUF 91) oe SHRP (Ey) 
COMMONS SHARPS, XIACe SHAPE (6520) 05G (3039) oSW (493) 0K (3020) LD 1120) 
EGULVeLENCE (SRP, SHAFE) 
NENFL ONC TM 
Co 2hu ITiglINnt 
SS#stuwlyer)) 
TT=sTuw leet) 
ULESTUW (ae ]) 
wlesTumlael |) 
CALL HRICRKR Ee (SSeTT UU sh IM,NEL INSIDE) 
DVOL Sas K GAC 
FURCEA SAC OW, 
AS20 
PO SU JS loNeEL 
CCBSHAPE (495) OR IAC 
00 SO 1®leNUIM 
WSense) 
KP O(N?) &Sh ape (19d) 
FORCE (Te J) SF ORCE (i 94) CC 
NS 70 
NO 15: Jwelynr 
CCBSrtr (4) fOvUL 
Og 1uy I@ley 
STcpenSpe3lepenS) «cc @SHP (yy 
NSBENSeY 
CONTINUE 
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Lh1460C 
Mrale 
FimyAze 
FLY AIC 
FLM144C 
FLMIA9C 
FLM)AOC 
ELM14A7C 
Fy ete 
FLRIAIC 
FLY) 79C 
Fim) ?lc 
FLM) 72C 
FLY) 73C 
FLMI74C 
FLM),75C 
FL™) 76C 
FLMI77C 
FLm178C 
FLM}79C 
FL™)20C 
FLY) RIC 
Fuwy?ec 
FLK183C 
ELM) BGC 
ELM} aSC 
FLY) 86C 
FLMITR7C 
FimyPeC 
wot fe 
VSI 2¢ 
vs! 3C 
vSt 4¢ 
vST Sc 
vst oC 
vst 7¢ 
VST 8¢ 
vST Cc 
vst 10c 
VST 11C€ 
vST 12C 
vST 43C¢ 
VST 14C 
vST 15¢ 
vST 16¢ 
VST 17e 
vSt 16C 
9c 

vel 
vST alc 
vst 2eC 
Yost 23¢ 
VST 2ec 


vs 

vS 260c 
VST 7c 
vst 28¢ 


Coes 


luo 
eoa 
elu 
220 
230 


400 


419 


tol 


tol 


tue 
103 


CLE PTOI GE AITO ot 5 


RE Th, 
ENL 


SLRROUTINE aSSEMBUNLIM a NDE eNEL oNSTFE SST SE STIF eFORCE 6 GGe KX eRO) 
Coecee oogtmee POOOHH 7/1) 3/7 I SHOHHHHESEEHHDEPHESHEREEHHHEEHUNEOOERSE 


DIMENSTON ESTIFCNSTH ot STF) SE C1) eFORCE C1) 
BSSEMELY Of ELEMENT INTO THE STIFFNESS 

RL @ REL ONOIM 

NS = (NUM CNLOLIIZE 

J) = appre (Neled) 


NDL_ & NUIM ¢« } 

po 23, y * LoNEL 
HO 22: N ® jye_NiIM 
2) = WO TON 

1] = yl 

bo 21 I ® yoNEL 
ay 


IF (Ter Qey) GC EN 

CO 2Uy ™ ® LeNUIM 

Ml 2 anl * w 

CC = dTI(Nd) 

ST{NS) F 900 

FSTiRCTPeM; pur eNp) = ESTIFCILOM) suieNd) © XFecc 
cc = GG*CCc 

Ect IF ilvenpesremy) & EstIFCLLenteuyemy) « CC 
TF (MpeoNbeN)) GO TU eUD 

DU you ® = yeNDIM 

FerIF il iemelp eh) S ESTIF CTL eRest en) « CC 
NS SWS: 3 

Th i@ PY = NYE 

CONT IAVE 

JIo= yl = NUF 

WO S wELON GS 

00 30u Jb * eahu 

Ye? me 

DO 360, 1 ® Jow 

ESTIF (Teg) = FSTIRF (Ley) © ESTIF (uel) 
ESTIF(JeT) = ESTIF EI eu) 

00 40u Telepnu 

FoTIF cdot = ESTIF (lol) exO@FORCE CL 

KE TUR 

END 

SLBHOLT INE [NSHAr(SeSGeShaPE oNDIM) 
DIMEN EON SFAPE (aed) es SO69-¢3) 0563) 
bO 16) L * penOIM 

Share (4e1) = le. = S(])##2 

DO yoy J ® pedDIM 

SHAPE (To J) & MeeC@SUtusleSly) 

KE TUR 

FRO 

SUBROUTINE [NCOMF (Ast eCoSoF oME sNE OT Sw) 
CIMFINSTUN ay DeMEy eB ( GedhE) oC (MED eS (NEG NE) oF (NE) 
GC TU fau8200910u)» TS* 

KEOUCE ARKarS 

OU 105 nelowe 

pe * Bingen? 

If (Peep ueren) GO Tu 105 

IF (NerWemME) GO TU lus 

nee ed 

HO 1¢¢ Leth Me 

BP = wtleN)/P 

TF (AR Eee) GO 10 103 

Ctl) & Cll)y © AP@CIN) 

DO yo, JelewE 

Ailes) = Ailes) © Aimed) @Ap 

Aljel) * Alpoy) 

00 Jue JELeoneé 

Bitegs = Pipeur = Bineur ep 

CONT Te 

IF (1SweNFely GO Tu 1U5 

OG faa EF # pont 

AP © glivgli/P 

F(l) = FOL) = APeC(W) 
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3)1 


StieJd) = Stf{sy) = 


rp 
« 

‘ 
uv 


= leweé 

= yq9nl 
Chak, <F: 

Lic tt 


V)sAINON 
RETUHIN 


[eve 


0033) 
tsTiF ¢ 
J>oeuks 


Clistte 


Lot INT 


2 WEL « KOM 
J = |oNELP 


< 
atCuwurrw 


\ 
prAPE LC) eu) #D 
SHAPE (e095) *D 
Chy Oxy 
Cheeks 
CBN 
Cyaan 
CIPO 
Clle@yn 


> 

~~ 

nx 
enwenw ot 


IF (NU |More?) 
fh SHePE (GQ ud oO 
Al3 ® Cle*Zén 
A4> = CRIM 
ASI = CliV4en 
60 303 | = JeNDIM 
Anh = 


GO 


mMPATIb6LE 


) 


HINe J) PAP 


DISPLACEMENTS 


(1,9) eF (U) 


y * AUNs J) ®C (J) 


LNSTIF CEST IF eFORCE sSGUeSTUMSLINT 9 NDF eNDIMONELONSIDVE 9 MMe 


Gol) sFORCE (1) »SGU (993) oSTUM (4,1) 


v 
v 


| 
v 


QO 301 


SHAPE Cy oT @NEL) 


YR 2 gMAPE (oe IT *NEL) 


IF OND L4e Fee 3) ZN 
EstyFati vol ) 
ESTIF CIR eyret? 
FSTIF (LT lels yt 
EctilFilieteurel) 
IF (NDI MeF lee) GU 
ESTIF (IE eyleed 
FsTIF(ti*ieyite) 
FSTIF(lieesst ) 
FSTIF CL bees ytel) 
ForlF(lie2eyie2) 
I] = 4) * WUF 

J) @ gt # "UF 
CONTI Nut 

RE TUR 

FRU 


SUBRUUTI NE PSTRES(NUIMe SIGsP 1 sP29P 3) 


DIMENSTUN S[G0%) 


SHAPE (3e1¢NEL) 
EstiF(ry ol ) 
ESTEECE? eS *P) 
ESTIECLleleyl ) 
Cor lF(lyelesied) 
O 303 
ESTIF (TI) syl*e) 
EsrIF (1itleylt2) 
ESrTlF(lle2e4l ) 
ESTIF (rleesylel) 
Esrik (lieesulee) 


+f oe 


eer ee 


AN@AL] 
AN@ale 
xN@®ael 
xN®A22 


AN@® al] 
Yn®413 
£N@A31 
YN@a42 
YN*®A2) 


e+ oe 


coe ee 


IAC eSPaAFF (4920) 9561363) oSKR0 393) oh (9020) ol O(120) 


HIKE (STUW (LoL) 6 STUm (29) eSTUW(39L) pNOTMyNEL sNSTOE) 
LNSr ah (STUN (ghd pSG0g SHAPE (1 NEL %]1) »NDIM) 
alUwlaelL ISSAC 


Yn@4o) + dnedao 
YNe@a22 
YN@A31 
YN®A32 © Znedae 
ZN®A22!) 
ZnN*®A21 
XN@AG2 
Zne@ale 
ZN®453 © KN#A22 


DATA BE27 1 6416213956257 3097 9P12372 00943951 02393217 


STHESSES MUST we STURFD IN akRaY S1G6(6) 
| RUMAK TAUEKY TAU aKL yg TAUST Ye TAUHVesTAURL2 FOR 3° PROBLEMS 


TF (NDIM*L Tea eOR eNUIMOGT 3) 


TF CNULMeb ue 3) GU 


RETURN 
0 100 
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In TRE ORDER 


INC 
INC 


INC 
TNC 
INC 
NC 
Inc 
INC 
INC 
INC 
INC 
INC 
TNC 
TNC 
INC 
INC 
INC 
INC 
INS 
INS 
Ivs 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
rvs 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
INS 
pst 
eS 
PS 
pst 
pest 
est 
PST 


——— 


Cores 


Creve 


euu 


COMPUTE PRINCTRLE STRESSES FUH 2° PROKHLEMS 

Teme = 2265/ATAN(1,9) 

XIL # (SIGE) * S1G43))42, 

Rl2 = (STGty) © SIGL3)IF26 

HRO = SURE (xT ePATe * SIG(2)*SIG(e)) 

Pl] = ATL © RRO 

P2 2 AT) & BHO 

PI = srel 

TE (Xd ee Nf 09) PS & TE APRATAN2(516(2) eAl2) 

RE TUR: 

COMPUTE FRINKCTIPLE STRESSES FUR Jen PROZLEMS. 

RrU S De 

Xi) = C31Giy) © FiGla) © STG(O))7360 

Klee = SIGE) ) *CSITG 14) *S1G616)) *S1G6 (4) #51616) S162) ##2-51G (3) ##2 
i =SIG6(5) **2 

ALTQ = SIGC YY PSTG 14) F916 (6) 22451612) S16 03) #91615) 
t *SLG CL) #81609) #426816 (4) 8516 (3) 082251616) S162) o#2 
DI2 = 3.ex];exl) - x12 

TauUG * us 

TF (D1e ete%,9) GO TU 200 

Vauv = Sgr (201 e736) 

Teme = (XL 3 6 (OTe © KEL@AP LI OAT PRIZs TAYE e3 
BHO = SURT (CaaS (hLesTEMPaTEMP) ) 

RRL #® ATAN2 (HOS TEMB) 736 

COMPYTE any RETURN PRINCIPAL STRESSES 

Pl = ATL © PAUIPRI2"®COS (RHO) 

be = ATL * fauUIPRI2*COS (RHO * PIe3) 

P3 = APS * FAUDPRIAPCOS (RHO * PI23) 

RETURN 

ERD 
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pst 
Pst 
pst 
est 
PST 
Pst 
est 
pst 
Pst 
PST 
PST 
pst 
eS! 
Pst 
eST 
est 
est 
PST 
PST 
est 
Pst 
est 
pst 
psi 
est 
pot 
est 
est 
PSt 


KC 
9¢ 
joc 
11¢ 
12¢ 
13C 
14¢ 
19¢ 
19¢ 
L7C 
18§ 
c 
20C 
alc 
2ec 
3c 
24 
25 
26C 
2qc 
abc 
29C 
30¢ 
31¢ 
32C 
33¢ 
34C 
35C 
36C 
37C 


